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ABSTRACT 


High Energy Lasers (HELs) are a new class of weapons that may be of great 
value to the Navy in the near future. A high-power Eree Electron Laser (EEL) is being 
designed using short Rayleigh-length resonators to increase the spot size at the mirrors 
and hence avoid mirror damage. Three-dimensional simulations are used to study the 
effects of an electron beam misalignment (electron beam tilt). This thesis shows that the 
proposed design is tolerant of typical electron beam misalignments. The performarce of a 
step-tapered undulator is also studied for the 100 kW proposed upgrade of the Jefferson 
Laboratory EEL. The results of this research show that the gain is above the required 
threshold for the 100 kW design while the energy spread does not change significantly 
over any undulator design. The spectrum of the proposed EEL shows that most of the 
power is concentrated around the fundamental frequency. It is shown in this thesis that 
smooth EEL pulses can significantly reduce the negative effects of absorption and 
scattering. Recent HEL science and technology developments are discussed for both Eree 
Electron and Solid State Lasers. 
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EXECUTIVE SUMMARY 


High Energy Lasers (HELs) have two characteristics that make them particularly 
valuable for ship self-defense: they are extremely fast and extremely precise. In addition 
to the potential of using high-energy lasers against moderately hard targets, laser systems 
could have the potential for lethal precision engagement against classes of soft targets, 
including personnel and light vehicles. The high power Eree-Electron Lasers (EELs) and 
Solid-State Lasers (SSLs) are the most likely candidates for developing a Navy high- 
energy weapon system. 

SSL technology, including fiber-laser technology, is often presented as the 
enabling technology for new HEL tactical missions. However, thermal control is the bane 
of all SSL systems, since the increase in temperature distorts the lasing medium. Thus, 
efficient thermal control of laser media needs to be pursued. Of all the solid-state laser 
concepts, probably the closest to achieving 100 kW level output power is the heat 
capacity laser (SSHCL), where lasing and cooling occur separately to avoid the growth of 
temperature gradients across the active area preserving beam quality. 

EELs can produce powerful, highly coherent, optical beams but also place 
extreme demands on the resonator mirrors. A critical challenge for EELs is to maintain or 
even improve electron beam brightness as the beam current is increased. Another concern 
for EELs is making them compact and rugged enough to meet shipboard packaging 
requirements. 

A high power free electron laser for ship defense is being designed at the Naval 
Postgraduate School in collaboration with Jefferson Laboratory. A short Rayleigh length 
resonator has been proposed to increase the spot size at the mirrors, and hence reduce 
mirror damage while meeting shipboard packaging requirements. This implies a nearly 
concentric cavity with a very small optical mode waist. A slight misalignment of the 
electron beam, either an offset or a tilt, could conceivably reduce the overlap between the 
electrons and the optical mode in the center of the undulator, resulting in less gain and 
efficiency. Three-dimensional simulations showed that the proposed design is tolerant of 



typical electron beam misalignments (either an electron beam tilt or an offset), providing 
enough efficiency for the MW level weapon system requirements. 

One transitional step toward the 1 MW output power required for a laser weapon 
is the proposed 100 kW upgrade of the Jefferson Laboratory PEL. Additional simulations 
were conducted to explore the use of a step-tapered undulator, which alters the resonance 
condition halfway through the undulator, in order to improve the EEL performance. The 
results of this research showed that step-taper is not always so effective; towever, in all 
cases explored, the weak field gain was found to be above the required threshold of 21% 
for the 100 kW design. The induced energy spread at saturation in strong fields did not 
change significantly over any undulator design. 

The spectrum and the power distribution of the proposed EEE output were also 
studied. Approximately 90% of the total power of the laser beam is concentrated in the 
main lobe (around 1 pm wavelength), while significant amount of power is contained to 

higher and lower harmonics due to short pulses. Atmospheric absorption and scattering 
were also studied for a specific maritime aerosol model showing the amount of power 
that is delivered through the atmosphere over a distance of 5 km. 



I. INTRODUCTION 


The primary Navy mission for a modern weapon system is cruise missile self 
defense. Cruise missiles have increasingly sophisticated capabilities and represent a very 
significant threat to present and future Naval operations. The supersonic cruise missile 
threat to Naval vessels, such as the SS-NX-27 Novator Alfa is well known and shown in 
Figure 1.1. The potential engagement time is limited given the distance to the horizon and 
speeds greatly exceeding Mach 1. 

This threat is compounded by the post-Cold War mission shift of Naval 
operations from predominantly “blue water” activities to now include littoral 
environments. The greatest evidence of this shift is the “Littoral Combat Ship” (LCS) 
concept, which has been adopted and is under development by the U.S. Navy. This 
project involves the building of a new type of ship with warfighting capabilities 
optimized for versatility in the littorals. It must have self-defense systems enabling it to 
survive threats such as small boats, patrol craft, low-slow flyers and shore artillery. The 
priority of this mission is enhanced by recent terrorist actions, such as the USS COLE 
incident, and ongoing terrorist threats. A high-energy laser (HEL) weapon system is an 
excellent candidate to fulfill this diverse role. 



Range 


220 km 


Speed 


Mach 0.6 to 0.8 subsonic 
Mach 2.9 supersonic 


Payload 


200 kg 


Guidance 


GPS, Inertial and active 
radar homing 


(UNCLASSIHB)) 


Length x 
Diameter x 
Wingspan 


8.22 m X 
0.533 m X 
3.1 m 


Eigure 1.1. SS-NX-27 Novator Alfa. (Erom [Ref. 1]) 


Two characteristics of the HEEs make them particularly valuable for ship self- 
defense: they are extremely fast and extremely precise. The laser begins its attack within 
seconds of detecting its target and completes its destruction a few seconds later. In 
addition to using high-energy lasers against moderately hard targets, laser systems could 
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have the potential for lethal preeision engagement against soft targets, such as personnel, 
light vehicles, communications, etc. Overall, HELs have the potential to change future 
military operations in dramatic ways and add a new dimension to a wide range of 
missions. See Figure 1.2. 




Figure 1.2. Naval Threat Suite. 


Greece, the only state in the Balkans and the Eastern Mediterranean that is a 
participant in both NATO and the European Union, is situated in a very difficult 
geographic position and is surrounded by sources of tension and engagement. The 
Hellenic warships must be capable of operating under any possible environmental 
scenario as part of NATO forces, including open seas and coastal environments. 

Moreover, the peculiar geographical structure of Greece, composed of extensive 
coasts with numerous natural ports and bays as well as a great number of islands and 
islets, makes the coastal Naval operations absolutely essential to national security. 
Therefore, all of the previously mentioned threats are relevant to Greece. Thus, the idea 
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of using Directed Energy Weapons (DEW) appears very appealing to Greece and the 
Hellenic navy. 

At this point, some critical factors must be taken into consideration. The prime 
power requirement for such a weapon is estimated to be of the order of 10 MW, which is 
far beyond the power available from the Hellenic warships (3 MW). Also, the possible 
cost is expected to exceed $50 million [Ref. 2]. Especially at the beginning of production, 
this cost is much greater than conventional weapons. 

However, even though this may presently be unfeasible, the DEW concept may be 
a relevant and serious challenge for Greece and other countries during the next few 
decades. 

The high power Eree-Electron Easers (EEEs) and Solid-State Easers (SSEs) are 
the likely candidates for developing a Navy high-energy weapon system. The U.S. Navy 
is closely watching the continued efforts of the Army and the Air Eorce in chemical-laser 
source development, but intends to focus its source-development efforts on electric- 
powered laser devices, such as solid-state and free-electron lasers for the following three 
reasons [Ref. 3]: 

• the inherent dangers associated with the storage of chemicals on board 
ship 

• the need to develop a source at the appropriate wavelength for maritime 
propagation 

• the Navy’s recent decision, early in 2000, to adopt an electric drive 
propulsion system for DD-21 and follow-on ships 

Though both offer the benefits of using electrical pumping mechanisms, they each have 

challenges unique to their design. The solid-state lasers will need to overcome the 

thermal loading that a multi-megawatt system will place in the laser medium. The free- 

electron laser source will require high average current injector work, as well as work on 

high-power optical resonator mirrors and an electron beam transport system. 

There are also several challenges related to atmospheric effects and propagation 
in developing a high-energy laser weapon system that will operate in a maritime 
environment. The atmosphere just above the ocean is very different from the atmosphere 
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over land, and will require the selection of the appropriate wavelength, power level and 
pulse format in order to deliver a lethal amount of energy to the target. 

Absorption scattering and thermal blooming effects are the primary causes for 
loss of energy in this environment. See Figure 1.3. In reality, maritime conditions at a 
particular geographic location and at a particular time may favor different wavelengths on 
different days. This variable nature of maritime conditions provides a significant 
advantage to the free-electron laser. The ability of a free electron laser to select different 
wavelengths is an important capability for a maritime laser. 

Thermal Blooming Scattering 




Figure 1.3. Thermal Blooming (Due to Atmospheric Absorption) and Scattering 

Effects. (From [Ref. 1]) 

Chapter II of this thesis describes the current initiatives as well as science and 
technology needs to achieve weapon-class high-energy lasers and presents an overall 
comparison between the two technologies. 

Chapter III of the thesis gives an overview of the theoretical background and 
discusses the physics behind the Solid-State-Laser and the Free-Electron-Laser. 

Chapter IV describes the results of simulations for the 1 MW EEL. This work was 
presented at the 24th International EEL Conference in Chicago (9/2002) [Ref. 4]. 

Chapter V describes the results of simulations for the 100 kW upgrade Jefferson 
Laboratory EEL as a transitional step for the development of a 1 MW which were 
presented at the 23rd International EEL Conference in Darmstadt Germany (8/2001) [Ref. 
5]. 
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Chapter VI presents the power distribution of a short pulse PEL beam over the 
frequency spectrum and studies the effects of atmospheric absorption and scattering in 
the performance of EEL as a weapon system. 

Chapter VII summarizes the conclusions of this thesis. 
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n. HEL WEAPON TECHNOLOGY 


Solid-state-laser technology, including fiber-laser technology, is often presented 
as the enabling technology for new HEL missions, and in particular, tactical missions. 
However, thermal control is the bane of all SSL systems since the increase in temperature 
distorts the lasing medium. Thus, efficient thermal control of laser media needs to be 
pursued. Other technology efforts for high energy include the design and manufacturing 
of reliable diode pump lasers, and the phasing of laser modules. 

The Tree-Electron-Laser (EEL) is an all-electric laser that is a competitor to the 
solid-state laser. A recent demonstration of its 2.1 kW average power [Ref. 1] has put the 
EEL in the same category as the state-of the-art solid-state lasers. To obtain more power 
from an EEL, the energy in the electron beam only needs to be increased by increasing 
the average current. 

Theoretical calculations have shown that EELs have the potential to scale to 
multi-megawatts. There are, however, many technical challenges to developing a high 
power EEL. A critical challenge is maintaining, or even, improving electron beam 
brightness as current in the beam is increased. Another concern for EELs is making them 
compact and rugged enough to meet shipboard packaging requirements. 

A. FEE DESIGN 

A Eree Electron Laser (EEL) is a device that extracts kinetic energy from a 
relativistic electron beam and converts it into electromagnetic radiation. This is 
accomplished by sending relativistic free electrons (not bound to atoms) through a 
periodic magnetic field produced by a device called an ‘Imdulator” or “wiggler”. The 
spatially periodic magnetic field induces transverse forces on the electrons and causes 
them to produce electromagnetic radiation in the forward direction of the electron beam. 
Eree Electron Laser (EEL) oscillators couple out a small portion of the radiation with low 
gain per pass and store the rest. EEL amplifiers couple out most of the radiation with gain 
per pass. EELs have been demonstrated to provide an electrically driven, powerful source 
of wavelength selectable, coherent electromagnetic radiation. 
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Figure 2.1. Typical Elements of an EEL (Regenerative Amplifier). (Erom [Ref. 6]) 

The basic building blocks of a typical FEE as illustrated in Eigure 2.1 are 

• A source of free electrons (Photoinjector) 

• An RE accelerator (Einac), which accelerates the electrons to nearly the 
speed of light 

• A wiggler, using a periodic magnetic field so the electrons emit pulses of 
electromagnetic radiation 

• An optical resonator which outcouples a small portion of the radiation and 
recirculates the electron beam 

• An electron beam decelerator or dump, which safely captures the 
remaining energy from the electron beam 

Very significant progress in EEE technology has been achieved over the past five 
years and a credible path to a MW-class weapon system demonstration is being defined. 
The critical component issues for the EEE development are high-average-power injectors, 
high-current electron beam transport, high-power oscillator and amplifier designs, and 
high-peak-power density in the optical resonators. 

B. FEE CRITICAL COMPONENTS 

I. Injector 

The purpose of the injector is to produce short, periodic bursts (pulses) of 
electrons and feed them into the accelerator. During the pulses, typical peak injector 
currents are in the 1 to > 100 amp range with pulse widths ranging from 1 - 50 ps [Ref. 
7]. The average power of the electron beam is increased by increasing the charge in a 
pulse, the repetition rate of these pulses, or the energy of the electrons. The average 
current is given by 



















where q is the charge per pulse, To the pulse duration, T is the pulse period and Q. = \/T 
the repetition rate. To achieve megawatt level average power from an PEL, the injector 
current (charge per bunch x repetition rate) must be in the range of a half amp to a few 
amps. Thus, the development of high power injectors seems to be a key issue for a high 
power EEL. 

The injector energy is not recoverable in the accelerator. Erom this point of view, 
low-energy injectors are preferred. To date, most have been thinking in terms of a 10- 
MeV injector [Ref. 8]. Given an average current of 0.375 A, the required RE power is 
correspondingly 3.75 MW, which is a large amount of un-recoverable energy for 10% 
wall plug efficiency (~10 MW) in an EEE design. There is a payoff in principle by going 
down to, say, 5 MeV, in the injector. 

Due to these tradeoffs, hardly any dispute exists within the EEE community that 
the key development required to scale from the current kilowatt level to the baseline 
megawatt level is the high current injector. Additional research is required to ascertain 
the optimal energy of the injector. 

Injector options include a photocathode gun with DC acceleration, a photocathode 
gun with RE acceleration at room, cryogenic, or superconducting temperatures, or a CW 
thermionic injector with RE acceleration. A trade study is required to assess the relative 
merits of DC guns and RE guns for high-average-power IR EEE's against consistent 
standards, and thereby, identify the break points, such as bunch charge, between the 
technologies. In general, RE guns may have an advantage in terms of beam quality at a 
higher bunch charge [Ref. 6]. 

Presently, the Jfefferson Eaboratory kilowatt EEE uses a photocathode injector 
with a 350 kV DC (electrostatic) accelerator. It provides about 5 mA average current in 
60 - 110 pC bunches [Ref. 7]. The DC-gun performance is still probably sufficient for 
IREEE purposes up tD 1 nC, but it is roughly a factor 3 worse than could be achieved 
with a RE gun [Ref. 8]. 
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2. Accelerator 

In a linear accelerator, the bunches of electrons are accelerated to near the speed 
of light (relativistic velocities) using intense electric fields. The energy acquired by 
electrons in an accelerator is measured in electron volts (eV) which is the energy gained 
by an electron when it passes through an electric potential of 1 volt. 

When an electron is accelerated to near the speed of light, its kinetic energy is 

0 0 0 7 /O 

KE =(Y -1 )mc where the Lorentz factor is y = (l-D / c )' , t) is the electron speed, 

c is the speed of light in a vacuum, and me is its rest mass of about 0.5 million electron- 
volts (MeV). Typical accelerators for FELs operate in the range of 20 MeV (y ~ 40) to 
500 MeV (y -iOOO) [Ref. 7]. 

The most common electron accelerator for the PEL is a radio frequency linear 
accelerator (RE linac) which uses standing RE waves in the cavities to establish the 
accelerating electric fields. RE sources for EELs generally operate between 500 MHz and 
2 GHz at frequencies where high power tubes are commercially available. Electron 
bunches are injected into the accelerator at the RE frequency or at a subharmonic, such as 
every 10th cycle. 

The physical length of the accelerator is determined by the desired total 
acceleration (perhaps 100 MeV) and the achievable acceleration gradient, e.g., 6 
MeV/meter, of the particular accelerator. A goal of the accelerator community has 
always been to increase the accelerator gradient, and thereby allowing for more compact 
devices. 

Higher gradients are generally achievable using higher frequencies, while the size 
of the structure is also decreased. In addition, a higher RE frequency in the accelerating 
cavities is preferred because it enables one to keep the bunch charge smaller for a given 
average current and improve the beam quality [Ref. 9]. Roughly speaking, the beam 
degradation appears to be manageable for a sub-nC bunch charge, but becomes 
problematic at 1 nC where space charge effects degrade the beam unacceptably. 

However, the accelerating cavities store more energy and have reduced wall 
interactions with the electron beam at low frequencies. Lower frequency RE allows larger 
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cavities and improved engineering of water channels able to dissipate heat generated by 
ohmic losses from the RF energy more efficiently. This permits higher average currents 
and peak charges with less degradation of beam quality than can be achieved with the 
higher gradients associated with the high frequency cavities. The optimal frequency is 
unknown, but appears to be somewhere around 750 MHz for this MW design [Refs. 8 
and 10]. 

An alternative accelerator design uses superconducting RF technology that leads 
to a reduction of many orders of magnitude in resistive heating at the expense of adding 
cryogenic refrigeration. Jefferson Laboratory has demonstrated an FEL which produced 
2.1 kW of average power using a superconducting RF Linac to produce a continuous 
stream of micropulses [Ref. 8]. 

Superconducting radio frequency (SRF) accelerators with energy recovery are the 
appropriate approach in a system where the added complexity is offset by the need for 
long run-times at moderate to high average power and where system efficiency is 
important. Analysis of the specific mission application and candidate accelerator designs 
would be required to determine whether a room temperature or superconducting 
accelerator is the best choice. 

Over the past five years, SRF accelerators gradients have increased by more than 
a factor of two to ~ 10 MV/m, with 20 MV/m considered within near-term reach [Ref. 8]. 
Recent experiments have demonstrated greater than 1 kA peak current for high gain in 
amplifier systems with bunches of ~ 1 nC which are ~ 1 psec long. 

3. Wiggler 

The fee’s optical field is formed and amplified by the interaction of the 
relativistic electron beam with the spatially periodic magnetic field of the wiggler (Eigure 
2.2). Alternating the polarity of the magnets produces the desired magnetic field. This 
field’s strength can be adjusted by varying the gap between the fixed magnet poles or by 
varying the current through electromagnets. 

The wiggler’s period, , is the distance along the beam axis between magnet 
pairs. Wigglers are characterized by a dimensionless parameter K = (eBX^/2nmc^) - 
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9.345\ where B is the rms magnetic field in kilogauss over each period and \ is the 
wiggler period in centimeters. Wigglers are typically operated with ^ = 1 [Ref. 11]. 


Figure 2.2. FEL Critical Components. (From [Ref. 1]) 


The FEL’s output wavelength is given by X = (\/2y^ )(1 + ). Erom this 


equation, it can be seen that adjusting the wiggler period, the wiggler strength, or the 
electron beam’s energy can change the wavelength of the EEL. Eor operation at 
wavelengths in the near infrared (IR) or less, required beam energies are in excess of 10 
MeV [Ref. 9]. Eor instance, if = 1, = 3 cm, and the desired wavelength is X =1 pm, 


then the required energy is given by 


Y = 




170 ^ Eb earn — 85 MeV 


( 2 . 2 ) 


EELs also require electron beams with low emittance (small radius and angular 
spread) and very small energy spread (highly monoenergetic) [Ref. 7]. As the electron 
beam passes through the wiggler, only a small fraction of its kinetic energy is converted 
to optical energy. As energy is extracted, the electron beam bunches longitudinally and 
the energy spread increases. If the EEL uses energy recovery, only minor degradation of 
the electron beam is tolerable. Eor recirculating EELs no more than a few percent of the 
beam energy can be extracted and converted into optical energy each pass. If, however. 
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the remaining beam is to be disearded (as in a regenerative amplifier), then mueh more of 
the beam energy might be extraeted. This is usually aehieved by tapering the wiggler 
period or magnetie field strength to maintain the resonanee eondition as the eleetron 
beam loses its kinetic energy to the light wave. 

4. Resonator 

In an oscillator PEL, the optical cavity is oriented along the wiggler axis and 
extends beyond the wiggler length, bounded by the resonator optics as in a conventional 
laser (Figure 2.2). During successive passes of the optical pulses through the wiggler, the 
optical beam is amplified and a fraction of the coherent radiation is allowed to escape and 
become the outcoupled optical energy from the EEL. 

For maximum coupling efficiency, the optical beam within the Esonator must 
roughly match the electron beam size [Ref. 11]. This is quite small within the wiggler, 
ranging from 100 pm to a millimeter in diameter [Ref. 7]. As a result, the laser beam is 
very small and intense when it strikes the resonator optics. This problem is exacerbated 
by increasing the FEE’s average power. If resonator length is not a design constraint, this 
problem can be helped by separating the mirrors widely and allowing diffraction to 
spread the optical beam. For compact applications, i.e., to meet shipboard packaging 
requirements, this is not an option and alternate resonator concepts, such as a short 
Rayleigh length design or higher damage threshold optics, must be developed. 

C. MW-CLASS FEE PROPOSED DESIGN 

The FEE design that seems to be the most attractive for a MW class FEE weapon 
system is shown in Figure 2.3. The electron beam is created in a SRF photoinjector, 
which then introduces the beam into a recirculating loop consisting of an SRF accelerator 
and a cryo-cooled sapphire mirror oscillator connected by bending magnets. The high- 
energy-acceptance bends manipulate the beam phase space to preserve beam quality prior 
to reinsertion 180 degrees out-of-RE phase into the SRF accelerator. The electron beam 
is then decelerated and delivered to the beam dump at low energy. This key aspect of 
same-cell energy recovery, already convincingly demonstrated by the Jefferson 
Eaboratory IR FEE, leads to high wall-plug efficiency and minimizes the radiation and 
shielding requirements for the laser [Ref. 8]. MW-class FEEs are projected to be 10-20% 
wall-plug efficient. 
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Figure 2.3. Schematic Diagram of a MW Level FEL Weapon System. (From [Ref. 

12 ]) 

FEL technology has made great strides in the past five years with the recent 
demonstration of 2.1 kW at Jefferson Laboratory. FELs are now in the same regime as 
state-of-the-art solid-state lasers and, may in fact, be scalable to much higher powers. 

D. SSL DESIGN 

The typical structure of a Solid-State-Laser (SSL), which is similar to other 
conventional lasers, is shown in Ligure 2.4. The flash lamp pump inverts the electron 
population in the lasing material (rod), leading to energy storage in an upper atomic level. 
If this energy is released to the optical beam by stimulated emission, or the downward 
transition under the presence of light, amplification takes place. After spontaneous 
radiation (downward transition independent of the presence of light) is emitted along the 
axis of the laser, the system starts to oscillate if the feedback determined by the 
reflectivity of the mirrors is sufficiently large. The role of the resonator (optical cavity) is 
to maintain an optical mode whose losses are replenished by the amplifying medium 
through stimulated emission. 
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Figure 2.4. Elements of a Typical Solid State Laser. 


A considerable amount of progress has been made over the last 15 years in scaling 
solid-state lasers to higher powers. Kilowatt average power levels are now routine, and 
10 kW output, it is believed, will soon be demonstrated. Systems at the 100 kW level are 
being designed today, and conceptual designs for MW-level devices have been developed 
by multiple industrial sources. 

However, vastly increasing technology efforts are focusing on three keys to high 
energy [Ref. 6]: 

• design and manufacturing of reliable diode pump lasers 

• phased combining of laser modules 

• thermal control of laser media 

E. CHALLENGES FOR MW-CLASS SSL 

1. Laser Diodes 

Among the various light sources that have been employed over the years to pump 
solid-state lasers, the diode lasers are most promising for pumping high energy SSLs. 
The most basic form of a diode laser consists of a semiconductor in which one part is 
doped with electron donors to produce ntype material and the other part is doped with 
electron acceptors (holes) to produce p-type material. Application of forward bias in this 
p-n junction drives electron and holes into the junction region where electron-hole 
recombination takes place, releasing energy in the form of photons [Ref. 13]. 

Some of the major attributes of diode pumping are high pumping efficiency, 
improved beam quality due to the spectral match between the laser diode- laser emission 
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and the absorption band of the laser material, increased lifetime, directionality and a 
small emitting area of the diode output [Ref. 14]. 

However, laser diodes, as with all semiconductors, are devices very sensitive to 
temperature. Thus, the need to maintain the operating temperature within a relative 
narrow range requires a more elaborate thermal management system whe n compared to 
flashlamp pumped lasers. 

Diode arrays, arrays containing a number of stripe lasers on a single chip, have 
demonstrated more than 1 kW/cm^ output intensity [Ref. 6]. However, pumping 
HELs, designs which produce higher coherence, require high intensity pumping 
(> 10 kW/cm^) and thermal management is one of the principle challenges to achieving 

higher power diode arrays. Present diode arrays operate efficiently at junction 
temperatures up to 50-60°C. Operation at up to ~ 100°C has been demonstrated, but at 
decreased efficiency and lifetime. 

The peak power of the diodes also determines their packaging and cost. A typical 
engagement requires ~ 1 MJ of energy to the beam director, e.g., 400 kW output for 2.5 
sec. Since the laser medium itself will be only about 40% efficient, 2.5 MJ of diode 
energy must be delivered to the laser head. This corresponds to an average diode pump 
array power of ~ 1 MW delivered during the engagement burst. A pulsed mode operation 
of the diodes with 10% duty cycle would then require a peak power capability of 10 MW. 

Present cost is in the range of $10 per peak-Watt for pulsed arrays, $30 per W for 
CW arrays, and s dominated by the cost of the coolers. Lawrence Livermore National 
Laboratory (LLNL), Navy-Mantech, and Air Lorce Research Laboratory (ALRL) all have 
programs to reduce the cost and improve the reliability of diode pump arrays for high- 
energy lasers. 

2. Phasing of Laser Modules (Fiber Lasers) 

Another promising method to scale SSLs to high energy is to phase the laser 
modules. The “phase-locking” technology is an active research area in industrial, DOD, 
and DOE laboratories. Phase locking works towards coherently combining multiple laser 
modules. The Phase Integrated Laser Optics Technology (PILOT) program, conducted at 
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AFRL during the late 1980s and early 1990s, demonstrated coherent phasing of 900 low 
power diode elements with total output power ~ 10 Watts. 

Similar phasing techniques may be applied to current fiber lasers efforts to 
produce high-power fiber module devices. Current state-of-the-art fiber lasers produce 
diffraction-limited output but are limited in power to approximately 100 W [Ref. 6]. 

In principle, combining one thousand fibers, each operating at 100 W would 
produce a 100 kW laser system from an aperture less than 1 cm in area. The problem is 
that each individual fiber is an independent laser. As a result, the output of the individual 
fibers adds incoherently and the overall divergence of the fiber array is the same as for 
any one fiber. In other words, the array indeed produces 100 kW output power but in a 
highly divergent beam with a time-dependent, random coherence. 

If the technology to coherently couple the output of individual fibers, or the 
phase- locking of one fiber to another, can be developed then the promise of high power 
fiber lasers can be realized. The state-of-the-art laser is the achievement of two low 
power fibers successfully phase-locked together producing a coherent output of ~ 100 
mW [Ref. 6]. 

3. Thermal Control of Laser Media 

Thermal management is essential to solid-state laser operation since the increase 
in temperature distorts the lasing medium. The temperature gradients degrade beam 
quality and possibly deplete the laser gain by thermally exciting the lower atomic level. 
Both of these effects degrade the utility of solid-state lasers as weapons. Degradation of 
beam quality limits the achievable irradiance on target, while excitation of the lower 
atomic level limits the amount of laser energy and/or duration of the laser burst. Efficient 
heat removal usually dominates designs for high-average-power systems. 

The cooling techniques are described in Ref. 14. In most laser architectures, 
removal of the thermal load is accomplished by conduction through the laser material 
(convection on the surface) [Refs. 14 and 14]. The term continuously cooled lasers is 
used to refer to lasers where excitation and cooling of the laser medium occur 
simultaneously. These lasers have the longest history of achieving kilowatt-level output 
power. 
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However, the combination of volumetric heating of the laser material by the 
absorbed pump radiation and surface cooling required for heat extraction results in 
temperature gradients across the active area. This creates thermal lensing effects, due to 
a temperature-dependent variation of the index of refraction, which significantly degrades 
the laser’s beam quality. 

An additional issue associated with thermal loading is stress fracture of the laser 
material. Stress fracture occurs when the stress induced by temperature gradients in the 
laser material exceeds the tensile strength of the material. 

An incremental improvement to these effects can be achieved by using a “Slab 
laser with Zig-Zag optical path” (Figure 2.5). The rectangular-slab laser provides a larger 
cooling surface and essentially a one-dimensional temperature gradient across the 
thickness of the slab. Additionally in the zig-zag geometry, the optical beam does not 
travel parallel to the z axis. Instead the beam traverses the slab at an angle with respect to 
the x-z plane using total internal reflection from the slab y faces. Since the thermal profile 
is symmetrical relative to the center plane of the slab, the thermal stress averaged from 
one slab surface to the other (in the ideal case) is zero. However, despite the elegance of 
this design the performance of these systems has been noticeable poorer than predicted. 
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Figure 2.5. Schematic of a zig-zag slab Laser, (a) Slab geometry, (b) End view, (c) 
Optical path in y-z plane. (After [Ref. 14]) 


Diode-pumped Neodymium-doped Lasers (Nd:YAG, Nd:YLF, and Nd:YV04) 

originally dominated this class of lasers. Scaling to higher average power with 

continuously-cooled Nd-doped lasers may be difficult due to the relatively large quantum 

defect associated with Neodymium as explained in the next chapter. A possible approach 
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to achieving higher average power output from continuously cooled lasers is by using 
other materials, such as Ytterbium (Yb), in place of Nd, with lower quantum defect and 
require substantially less cooling and produce a better beam quality than a Nd:YAG 
device. 

F. HEAT CAPACITY LASER 

Of all the solid-state laser concepts, probably the closest to achieving 100 kW 
level output power is the heat capacity laser, which will be discussed in detail [Ref 6]. 
The heat capacity laser name is derived because lasing and cooling occur separately. See 
Figure 2.6. A pulse sequence begins by switching electrical current through the diode 
arrays, which excites the laser diodes and produces pump light for the laser crystals. As 
energy is absorbed within the crystals, the laser gain increases until it exceeds the cavity 
loss. At this point, lasing occurs within the cavity in free running mode. During the 
lasing cycle, pump radiation is continuously absorbed, while, simultaneously, heat is 
deposited in the laser crystal by the residual energy left behind. This energy is dominated 
by the quantum defect, which is equal to the difference in energy between the photons 
absorbed and those emitted ( ~ 0.36 eV for Nd lasers). 

The crystal temperature increases linearly throughout the lasing cycle. Laser 
emission is terminated by the user when target kill achieved, or when the lower atomic 
level becomes thermally populated, or when thermally-induced cavity losses exceed the 
available gain. At this point, current to the diode pump arrays is halted and a cool-down 
cycle begins. 

The amount of heat that can be deposited in the medium before lasing ceases is 
determined by the medium’s heat capacity. Since it is the total heat deposited which 
limits total operation time, heat capacity lasers are often referred to as having a “thermal 
magazine.” The thermal magazine may correspond to several target kills. Once the 
crystal has reached its heat capacity limit, it must be cooled before another cycle begins. 
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Figure 2.6. Heat Capacity Laser Concept. (From [Ref. 6]) 

The HCL laser module envisioned for the Army’s Enhanced Area Air Defense 
(EAAD) mission can operate for ~ 10 seconds before reaching its heat capacity limit. 
Since each engagement could last between 2 and 3 seconds, a total of ~4 target 
engagements could be achieved before the laser has to initiate a cool-down cycle. 
Depending upon the cooling technology, a cool-down cycle of 0.5 to 5 minutes would be 
required before the system is fully reset to room temperature. 

A state-of-the-art heat capacity laser has been recently completed that is 
flashlamp-pumped and reaches 10 kW. It has 500 J per pulse at 20 Hz pulse repetition 
frequency operated at EENE. See Eigure 2.7. While only 1.4 kW (140 J at 10 IE) 
average power has been demonstrated to date, the laser has recently exceeded the pulse 
energy design goal significantly by operating in a single-shot mode at over 1000 J per 
pulse. Replacing the flashlamp pump sources with diode lasers and replacing the 
Nd:glass laser disks with Neodymium-doped Gadolinium Gallium Garnet (Nd:GGG) 
crystals, is expected to improve a power-scaled version of this laser to achieve more than 
100 kW (500 J pulses at > 200 Hz repetition rate). 
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The High Energy Laser Systems Test Facility (HELSTF) received the 10 kW 
devices in August 2001 for integration into a testbed. Static tests are currently being 
conducted to evaluate and characterize the device in support of developmental testing that 
will lead to the 100 kW SSHCL prototype. 



Figure 2.7. LLNL 100 kW SSHCL Design. (From [Ref. 1]) 

There is a detailed Army/LLNL development plan for the 100 kW HCL 
demonstrator. See Figure 2.8. This plan is comprehensive in that it addresses most of the 
subsystems that will be required to achieve a viable heat capacity laser: diode arrays, 
laser crystal availability, and power supplies. Perhaps the most significant issue with 
heat capacity lasers that is not currently being addressed is cooling the active medium 
between lasing cycles. LLNL has suggested both gas (convective) and evaporative 
“mist” cooling. Little work has been done in this area, and it is unclear how these 
cooling methods will perform, including their impact on system complexity. However, 
the cool-down cycle is a critical factor that directly affects the reset time between 
magazines and may impose limits on the operational capabilities. 
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GA IS SYSTEMS INTEGRATOR & THERMAL MANAGEMENT SYSTEM 
DESIGNER FOR THE HEAT CAPACITY SOLID STATE LASER* TEAM 





Figure 2.8. SSHCL 100 kW Demonstration Concept for the Army. (From [Ref. 3]) 

G. CONCLUSIONS 

This chapter discussed the recent HEL science and technology developments for 
both FELs and SSEs. Areas needing improvements were addressed. 

In general, EEEs can produce powerful, high coherent, optical beams but also 
place extreme demands on the resonator mirrors while SSEs in general produce a larger 
beam spot and thus subjects optics to lower a damage threshold. They do not, however, 
achieve the desirable beam quality. Efficient heat removal from the laser media is 
needed in high-power SSE designs, while high-power EEEs require high-brightness 
injectors and novel optical cavity designs. However, considering the vastly increasing 
research development efforts, both technologies have the potential to meet the needs of a 
weapon-class laser. 

Eiber optic lasers and the Heat Capacity Easer project are examples of the 
readiness of solid-state laser technology for rapid maturation. On the other hand, the 
unique advantages of EEEs, such as wavelength selectability and tunability, high beam 
quality, the potential to achieve high peak power without distorting the lasing media and 
long run-time, could favor EEEs for weapon-class lasers in the near future. 
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m. PHYSICS OF SOLID-STATE AND FREE FEE CTRON LASER 


The word 'LASER' is an acronym standing for 'Light Amplification by Stimulated 
Emission of Radiation'. The output of a laser can be infrared, visible, or ultraviolet 
radiation in a pulsed or a continuous beam, with average power from less than a milliwatt 
to several megawatts of power. However, all conventional lasers have the following in 
common: 

• a lasing medium which can be a solid, liquid, or gas, 

• a means of pumping energy (optical, electrical, mechanical or chemical) 
into the lasing medium. 

In a Solid-State Laser (SSL), the active medium is a solid material and optical 
energy is used as an external pump source. The Tree Electron Laser (EEL) represents a 
radical alternative to all conventional lasers where the active medium is a relativistic 
electron beam. This chapter will outline the basic ideas underlying the operation of both 
types of lasers. 

A. PHYSICS BEHIND SOLID STATE LASERS 

To understand the operation of a SSL, which is actually similar to all conventional 
lasers, it is necessary to know some of the principles governing the interaction of 
radiation with matter. Many of the properties of lasers may be readily discussed in terms 
of the absorption and emission processes, which take place when an atomic system 
interacts with a radiation field. Einstein’s discovery of stimulated emission provided 
essentially all of the theory necessary to describe the physical principle of the laser. 

I. Einstein A and B Coefficients 

Consider an idealized material consisting of atoms with just two nondegenerate 
energy levels, Ei and E 2 , having populations of Ni and N 2 , respectively. Of particular 
interest is the interaction between an atom and a photon of frequency f = f 21 y such as 
hf 2 i = E 2 -Ej where his Planck’s constant. Three types of interactions can be identified: 
absorption, spontaneous emission and stimulated emission. 
a. Absorption 

An atom in state 1 absorbs a photon from the radiative field and thus 
converts the atom into state 2. The rate at which this process takes place depends on the 
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concentration of absorbing atoms and the intensity of the field from which they extract 
the energy, as is expressed by 


^ =+B,,p,(f)N„ 

V Jabs 


(3.1) 


where B 12 is a constant of proportionality and p^(/) is the radiation density. 
b. Spontaneous Emission 

After an atom has been raised to the upper level 2 by absorption, the 
population of the upper level 2 decays spontaneously to the lower level at a rate 
proportional to the upper level population, as expressed by 


dN, 

dt 


— A21N2, 


Jspont 


(3.2) 


where A 21 is a constant of proportionality and called the spontaneous transition 
probability (the probability that an atom in level 2 will spontaneously change to a lower 
level 1 within a unit of time). The reciprocal of transition probability of a process is 
called its lifetime, 

X2i=A-;. (3.3) 


Spontaneous emission is a statistical function of space and time, and there is no relation 
between the individual emission processes; the photons emitted incoherently. 

c. Stimulated Emission 

Emission takes place, not only spontaneously but also under stimulation 
by electromagnetic radiation of the same frequency. This process is the reverse of 
absorption. The atom gives up its excess energy hf to the field, adding coherently to the 
intensity. Thus, the added photon is at the same frequency, at the same phase, in the same 
sense of polarization, and propagates in the same direction as the photon that induced the 
atom to undergo this type of transition. The rate depends on the density of atoms to be 
stimulated and the intensity of the stimulating field. 


^ =-B2,pWn,, 

y Jstim 


(3.4) 


where ^21 again is a constant of proportionality. 


24 



If we combine absorption, spontaneous, and stimulation emission as 
expressed by (3.1), (3.2) and (3.4) we can write for the change of the upper and lower 
level populations 

()N ()N 

(3.5) 

At thermal equilibrium, the time rate of change must be zero. Therefore 

A,,N,+B,,pXf)N, = B,,pXf)N ,. (3.6) 


Using the Boltzmann equation and the Plank formula for the black body radiation law, 
the following relations between A’s and B’s which are known as Einstein’s relations are 
obtained [Ref. 16], 



and 

Si 


(3.7) 


where gj, g 2 are the degenaracies of levels 1 and 2 respectively, = c/n is the speed 

of light in the laser medium, n is the index of refraction and c is the speed of light in 
vacuum. 

2. Atomic Lineshapes 

In deriving Einstein’s coefficients, a monochromatic wave with frequency f 2 j 
acting on a two- level system with an infinety sharp energy gap hf 2 j has been assumed. 
The interaction between an atomic system having a finite transition linewidth Af and a 
signal with a bandwidth df will now be considered. The result is that a distribution of 
photon frequencies g{f) can be emitted (or absorbed). This relative distribution is called 
the lineshape function, with g(f)df being the probability that a given transition will 
result in an emission or absorption of a photon with energy between hf and h(f+df). 

The linewidth and lineshape of an atomic transition depends on the cause of line 
broadening. Optical frequency transitions in solids can be broadened by lifetime, dipolar 
or thermal broadening, or by random inhomogeneities. Two generic types of processes 
contribute to the width of a spectral line. 
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a. The Homogeneous Broadened Line 

The essential feature of a homogeneously broadened atomic transition is 
that every atom has the same atomic lineshape and frequency response, so that a signal 
applied to the transition has exactly the same effects on all atoms in the collection. 
Mechanisms that result in a homogeneously broadened line are lifetime broadening, 
collision broadening, dipolar broadening and thermal broadening. The lineshape of 
homogeneously broadening mechanisms leads to a Lorentzian lineshape for atomic 
responses [Ref. 17] as expressed by 

(3.8) 


gif) = 



(/-/o)" + 


2 

[2lt J 


V ^ y 



where is the center frequency, and A/' is the full- width half- maximum of the power 
spectrum. The factor Af/2K assures normalization of the area under the curve. 

b. The Inhomogeneous Broadened Line 

Mechanisms, which cause inhomogeneous broadening, tend to displace 
the center frequencies of individual atoms, thereby broadening the overall response of a 
collection without broadening the response of individual atoms. Different atoms have 
slightly different resonance frequencies at the same transition. As a result, the over-all 
response of the collection is broadened. Examples of inhomogeneous frequency-shifting 
mechanisms include Doppler broadening and broadening due to crystal inhomogeneities. 
The inhomogeneous-broadened linewidth can be represented by a Gaussian frequency 
distribution [Ref. 17] 


g{f) 


A/ 


ln2 


, 1/2 


V ^ y 


exp 


V J 


ln2 


(3.9) 


3. Amplification in a Laser 

The interaction between two linewidth-broadened energy levels with an energy 
separation centered at /„, and a half width of A/ characterized by g(f) and a narrow band 
signal with bandwidth df « A/ will be considered next. Equation (3.5) then becomes 


dt 




(3.10) 
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Using (3.7) 
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^ 21^2 
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Nx 
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(3.11) 


where = cpyn represents the intensity of the radiating field, is the radiation 
density of the field and a (/) = A 2 ig(f)'k^ IS!i¥ is the eross-section area of the beam. 

Sinee eaeh transition earries energy hf and ! ¥ represents the transition rate, 
we have 


-hf 


dN^ 

dt 


dz 


(3.12) 


where z is the direction the beam travels. 

Combining Equations (3.11) with (3.12) and neglecting the spontaneous emission 

we get 

|7(/)z 

If = I¥ ^ , ( 3 . 13 ) 


which can also be written as 




(3.14) 


where 



( g ^ 

N 2 -^N, 

= ^21 28 (f) 

( 8 ^ 
N 2 -^N, 


[ 8 x J 

Sun 

[ 8 x J 


is known as the laser gain equation. From Equation 3.15, it is obvious that, in order to 
have gain, it is necessary to have N 2 > (g 2 /g /) Nj and for = Si ^ ^2 > 

According to the Boltzmann distribution, in a collection of atoms at thermal 
equilibrium, there are always fewer atoms in a higher-energy level E 2 than in a lower 
level Ej. Therefore, the population difference Nj -N 2 is always positive and gain does 
not occur. The situation N 2 > Nj, which is required for a laser operation, is called 
population inversion and is never observed at thermal equilibrium. 
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4. Population Inversion 

This section is concerned with how the necessary population inversion for laser 
gain is obtained in solid-state lasers. Considerable understanding can be gleaned on how 
laser devices are pumped and how their population densities are inverted by studying 
some simplified but fairly realistic models. 

The discussion until now has been based on a hypothetical two-level system. 
The pumping and laser processes in real laser systems typically involve a very large 
number of energy levels. The main feature can be understood, however, through the 
familiar three-level or four level idealizations of Figures 3.1 and 3.2. 




Figure 3.1. Simplified Energy Diagram of a Three-Level Laser. (Lrom [Ref. 14]) 




Ligure 3.2. Simplified Energy Diagram of a Lour-Level Laser. (Lrom [Ref. 14]) 
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a. The Three-Level System 

In a three-level laser, most of the atoms of the laser material are initially in 
the lowest level 1 (ground state). Excitation is supplied to the medium by the radiation 
pump source producing absorption into the broad band 3. Most of the excited atoms are 
transferred by fast, radiationless transitions into the intermediate sharp level 2. In this 
process, the energy lost by the electron is transferred to the lattice. The difference 
between these two levels - E 2 is known as “quantum defect” and represents the 
laser’s minimum possible heat load. Finally, the electrons return to the ground level by 
the emission of a photon. It is this last transition that is responsible for the laser action. 

The basic disadvantage of a three-level system is that more than half of the 
atoms in the large-populated ground state must be raised to the metastable level E 2 in 
order to create inversion, which requires strong pumping. 

b. The Four-Level System 

The four-level system is characteristic of the rare earth ions in glass or 
crystalline host materials. The pump transition extends again from the ground state, now 
level Eg to a wide absorption band E^. As in the three-level system, the atoms will 
proceed rapidly to the sharply defined level E 2 . The laser transition, however, proceeds 
now to a fourth, terminal level Ej that is situated above the ground state Eg. From here, 
the atom undergoes a rapid nonradiative transition to the ground level. In a four-level 
system, an inversion of the 2 —> 1 transition can occur even with vanishingly small pump 
power so the high pump rate necessary in the three-level system is no longer needed. 

After this brief description of the energy levels, it is now obvious that the 
metastable level is of paramount importance for laser action to occur. The relatively long 
lifetime of the metastable level provides a mechanism by which the inverted population 
can be achieved. Most transitions of atoms show rapid nonradiative decay since the 
coupling of the internal atomic oscillations to the surrounding lattice is strong. In typical 
laser systems with energy levels, such as those shown in Figures 3.1 and 3.2, the 3^2 
and 1 —> 0 transition frequencies fall within the frequency range of the vibration spectrum 
of the host crystal lattice. Therefore, all these transitions can relax extremely rapidly by 
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direct nonradiative decay, emitting a phonon to the lattice vibrations, with lifetimes like 
X 32 ,Xio =10'^tol0" s. However, the larger 3—>0, 3—>1, 2^0, and 2^1 energy 
gaps often correspond to frequencies that are higher than the highest possible vibration 
frequency of the lattice. These transitions must then occur either by radiative (photon) 
emission or by multiple phonon processes. Since both these processes are relatively 
weak compared to direct single-phonon relaxation, the high frequency transitions will 
have much slower relaxation rates (Xjj ~ 10^^ to 10'^ s). Therefore, the various levels 
lumped into level 3 will all relax mostly into level 2. Level 2 is metastable (long-lived) 
since there are no other levels located close below it into which the atoms can be decayed 
directly. In order for the population inversion to increase, it is also required that 

"^30 32 21 10 • • 

To achieve the correct combination of relaxatbn times and the existence 
of broad pump bands, a suitable material must be used. 

5. Solid-State Laser Materials 

Materials for laser operation must possess sharp fluorescent lines, strong 
absorption bands, and a reasonable high quantum efficiency for the transition of interest. 
These characteristics are generally shown by solids (crystals or glass) which incorporate 
in small amounts of elements in which optical transitions can occur between states of 
inner, incomplete electron shells. The two principal elements of a solid-state laser 
medium are the host materials and the activator/sensitizer ions. 

a. Host Materials 

The host materials must have good optical, mechanical and thermal 
properties to withstand the severe operating conditions of practical lasers. Desirable 
properties include hardness, chemical inertness, the absence of internal strain and 
refractive index variations, and ease of fabrication. Glasses form an important class of 
host materials. The outstanding practical advantage compared to crystalline materials is 
the tremendous size capability for high-energy applications. The optical quality can be 
excellent, and beam angles approaching the diffraction limit can be achieved. Glass is 
easily fabricated and takes a good optical finish. Generally, the advantages of crystalline 
laser hosts over glass are their higher thermal conductivity, narrower fluorescence 
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linewidths and, in some cases, greater hardness. Some of the most useful laser hosts are 
the synthetic garnets: yttrium aluminum, Y 3 AL 5 O 12 (YAG), and gadolinium gallium 
garnet, Gd 3 Ga 50 i 2 (GGG). These garnets are hard, stable, optically isotropic, and have 
good thermal conductivities, which permits laser operation at high average power levels 
(SSHCL). 

b. Active Ions 

Among many others, the rare earth ions are natural candidates to serve as 
active ions in solid-state laser materials because they exhibit a wealth of sharp fluorescent 
transitions representing almost every region of the visible and near-infrared spectrum. 
Neodymium (Nd^^) was the first of the trivalent rare earth ions to be used in a laser, and 
it remains by far the most important element in this group. Stimulated emission has been 
obtained with this ion incorporated in at least 100 different host materials and a higher 
power level has been obtained from Nd lasers than from any other four-level material. 
Other rare earth ions in which laser action has been reported are Erbium, Holmium, 
Thulium, Gadolinium, Europium and Ytterbium. Diode-pumped Yb: YAG has especially 
become an important laser because it produces, due to small quantum defect, the smallest 
amount of crystal heating compared to any other major laser system. 

B. PHYSICS BEHIND FREE ELECTRON LASERS 

The Free Electron Easer (FEE) is a different kind of laser, which uses a 
relativistic electron beam as the active medium moving along a periodic magnetic field. 
The electrons are not bound to atoms, but they are not truly free since the wiggler field 
governs their motion. The FEE requires that the electrons undergo periodic oscillations 
in the transverse plane while moving at nearly the velocity of light along the z direction. 
The periodic oscillation implies an acceleration, which causes the electrons to radiate 
(spontaneously) perpendicular to their wiggling motion (like an oscillating dipole). If the 
z component of velocity is large ~ c (a highly relativistic beam), the radiation pattern 
is primarily in the +z direction and can be found from the relativistic Earmor formula 
[Ref. 11]. The spontaneous emission is stored in an optical cavity and amplified in 
succeeding passes by the relativistic electron beam (Figure 2.2). In order to understand 
how a relativistic electron can be coupled to the optical power, it is necessary to study the 
electron dynamics in the undulator. 
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1. Electron Motion 

The easiest way to derive the electron trajectories is for the case of a relativistic 
beam interacting with a helical undulator, or a circular polarized array of magnets. Near 
the undulator axis, the magnetic field can be expressed as 

= 5(cos(/c„z),sin(/c„z),0), (3.16) 


interacting with a circular polarized optical electromagnetic field given in cgs units by 

= ^(cosCvi/ ),-sin(\|/ ),0), (3.17) 

=£'(sin(\|/),cos(\|/),0), (3.18) 


where B is the magnetic field of the undulator, E is optical field amplitude, k^=2K/X^ is 
the wave number of the undulator, \|/ =(kz-Git -^) is the phase of the optical wave, and 
(j) is the initial optical phase at t=0, z=0. 


The problem being faced is to solve the relativistic Lorentz force equations 

d(ymc|3) 


dt 
d{ymc) 
dt 


-e(E+pxfi), 


-e^E, 


(3.19) 

(3.20) 


where the Lorentz factor is 

1 

Y= / ^ ^ 

Vi-PP 


(3.21) 


P -X) ! c is the dimensionless electron velocity, c is the speed of light, m is the rest mass 
of electron, e the electron charge magnitude and B = B^ + B^. 

Integrating Equation (3.19) using the undulator field from Equation (3.16) and 
the optical field from Equations (3.17), (3.18), and ignoring terms proportional to 1- P^ 
for relativistic electrons, we get 

p ^ = —(cos(^„z),sin(^„z),0), (3.22) 

Y 

where K — eB I kjnc^ is the dimensionless undulator parameter. 
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Equation (3.22) describes the transverse components of the electron trajectories as 
they pass through the undulator. It is now possible to understand how energy is 
transferred between the optical field and the electron beam. 

2. Energy Exchange 

The second Lorentz Equation, Equation (3.20), can also be written as 

= -gp .£ ^ , p j •(Ecos(v|/), -E sin(v|/ ),0) (3.23) 

at at me 


or 


— = (cos(fc^z)cos\i/ - sin(fc^z)sin\|/) 

dt mej 

dy eKE 

— =-cos(fc^z+tl/) , 

dt mey 


dy _ eKE 
dt mej 


cos((k +kg )z -tot +(|)). 


(3.24) 

(3.25) 

(3.26) 


It is convenient to define the electron’s position in terms of the electron phase ^ with 
respect to the optical field as 

^=(k+k^)z-(£>t, (3.27) 


so that Equation (3.26) then becomes 


dy 

dt 


= ^cos(C +(^). 
mey 


(3.28) 


The electron phase ^ represents the microscopic electron position within the 

electron pulse on optical wavelength, and is an important factor describing the energy 
exchange between the electron beam and the optical field. When the phase ^ +(|) of an 

electron falls between + 7t/2 , the cosine term will be positive, resulting in energy transfer 
from the optical field to the electron beam. The optimum energy exchange rate is 
described by a resonance condition. 

3. Resonance Condition 

As mentioned previously, the electrons entering the undulator are moving very 
close to the speed of light. However, special relativity requires that the finite energy 
electrons travel slower than the light beam. Resonance is described as the condition 


33 



where, over the length of one period of the undulator, one wavelength of light passes over 
an eleetron (Figure 3.3). 




At 




v/^ v OP 




< -> 


Figure 3.3. Resonance Condition. One Wavelength of Light (Traveling at Speed c), 
Passes over an Electron (Traveling at Speed ), as the Electron Travels Through One 

Undulator Wavelength. 


In a time interval At , a single electron with speed , travels distance \ = p^At 
and a co-propagating photon distance = cAt , so that 


^1-p ^ 


V. 


Equation (3.21) reads 


2 1 
y 02 0 2 


and making use of (Equation 3.22): 

y"(l+ii:") = l-p^". 

Assuming y » 1 and using binomial expansion yields 

1+K' 


Pyl- 


2y^ 


Combining Equations (3.29) and (3.32), 
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(3.31) 


(3.32) 


(3.33) 
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which demonstrates how the optical wavelength can be controlled by three parameters: 
the undulator wavelength, the undulator parameter (i.e., magnetic field strength) and the 
energy of the electron beam. 


In the PEL electron beam, there are many electrons uniformly distributed in all 
possible electron phases ^ as they interact with each wavelength of light. According to 
Equation (3.28), for a uniformly distributed electron beam, every electron, which 
contributes energy to the optical field, is balanced by another electron, which absorbs an 
equal amount of energy from the field. In order to understand this better, another 
important parameter, the electron phase velocity, must be discussed. 

4. FEE Pendulum Equation 

The rate of change of the electron phase is given by 

^ = (A: + A:Jp^c-a). (3.34) 

dt 


Taking another derivative and using Equation (3.31) for [3^,, we get 

dX 


dt^ 


■■{k+k^)c 


( 1 + kM 


y^ 

V ' 7 

[yj 


(3.35) 


Eor relativistic electrons, k» k^. Substitute y from Equation (3.28) and using the 
resonance condition from Equation (3.33), we get 


dt^ 


2eKEk„ 

y 


cos(^ +(|)). 


(3.36) 


Equation (3.36) can also be written in a dimensionless form 


V = 



=1 a lcos(^ +(|)), 


(3.37) 


where 


r o 


l = d(...)/dx, \a \=4'KNeKLE/y^mc^, a = \a\e‘^ 


is the complex dimensionless 


V J 

optical field, x=ct/L is the dimensionless time and = L[(k + k„)(3^ -k]is the 
dimensionless phase velocity (recall L = is the undulator length). 
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Equation (3.31) is called the FEL pendulum equation and describes an important 
EEL physical phenomenon called “eleetron bunching”. Whenever the total phase of an 
electron falls between ±71 / 2, the cosine term will be positive resulting in an acceleration 
of the electron. When the phase is between 71/2 and 37t / 2, the electrons decelerate. 
This change in phase velocity of the electrons causes the electrons to bunch. Under 
different conditions, just above or below Ksonance, the bunching occurs at different 
values of ^ . It is actually this relative location of the electron bunching that determines 
the net energy transfer to the optical beam. Eor maximum energy transfer from the beam 
to the optical field, “bunching” around the phase ^ +(}) =7t is required (according to 
Equation 3.28). 


It can be shown that the energy change can be related to the ehange in phase 
velocity as [Ref. 11] 


Av = 47C 


^ Ay ^ 
Y 


(3.38) 


5. Phase Space Analysis 

The microscopic motion of the electrons over an optical wavelength (bunching) is 
described by phase-space evolution of electrons. In order to understand how the electrons 
evolve in phase space, it is useful to simulate their motion over an optical wavelength 
based on the principal equation of motion, or the Pendulum Equation (3.37). Twenty 
sample electrons are initially distributed uniformly in phase from ^ = -7t /2 to 371 / 2 as 
they enter the undulator where the interaction occurs. 


In the beginning, initial eleetron phase velocity = 0 is used, which corresponds 
to the resonance eondition. In Eigure (3.4), it can be seen that half of the eleetrons move 
up in the phase spaee plot and gain energy from the optical field, according to Equation 
(3.38), and half move down and lose energy to the optical field. At the end of the 
undulator, the electrons are bunched around ^ +(|) =7t /2 where (j) = 0. The net gain is 

(G = O) and the optical phase shift is small. 
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Figure 3.4. Electron Evolution in Phase Space, Optical Eield Gain G(x) and Phase 
(p(x) for a Single Pass through the Undulator with Electrons Injected at Resonance 

(Vo=0). 

In Eigure (3.5), the evolution of twenty electrons with initial velocity = 2.6 
(above resonance) is depicted. In this case, the electrons on the left move faster than 
those on the right resulting in an electron bunching around ^+(t)=7t((t)~0) while they 

also fall to lower values of V , implying a net energy transfer to the light beam. The 
electrons still follow pendulum phase space paths, but significant energy extraction and 
gain now occurs, as seen at the right. 

The separatrix curve in both figures is drawn to separate the closed and opened 
orbits. Close orbit paths in phase space remain between -7t/2and3Jt/2, while open 
orbit path extend beyond these limits. These phase space simulations use periodic 
boundary conditions at -;t / 2 to 3Jt / 2. 
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Figure 3.5. Electron Evolution in Phase Space, Optical Eield Gain G(x) and Phase 
(p(x) for a Single Pass through the Undulator with Electrons Injected Slightly Off 

Resonance = 2.6 . 


After describing the electron’s evolution, it is necessary to derive the wave 
equation that drives the optical field amplitude and phase evolution. 

6. FEE Wave Equation 

The optical wave equation can be derived by starting with Maxwell’s wave 
equation in the Coulomb gauge [Ref. 18] 
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2 1 ^ 

V- 

c' de 


2 A 


4 471 - 

A =- J L, 

c 


(3.39) 


where A is the vector potential and /lis the transverse electron beam current. Eor an 
electron beam interacting with circular polarized light in a helical undulator, the optical 
vector potential is [Ref. 11]: 

(sin\|/,cosV|/,0), (3.40) 

k 


and the electron beam current is 

7;;=-ecM'(x-;;), (3.4i) 


where 6^(x)is the three dimensional delta function and r. is the position of the 4i 
electron with charge e magnitude and transverse electron velocity [3 given by Equation 
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(3.22). It can be further assumed that the optical field amplitude and phase vary slowly in 
time and space as that 


^«kE. 

dz 

dE 

— «(OE, 
dt 

(3.42) 
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dz 
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dz 

(3.43) 


Then, Equation (3.39) can be written 

(dE IdE^ fdd 194)3 47t- 

2 —H-— (cosv[t,-sinv|/,0)-2E ——i-— (smv|/,cosv[t,0) =- J± (3.44) 

^oz c at J yoz c at J c 

Use the two unit vectors Ej =(cos\|/,-sin\|/,0) and Ej =(sin\|/,cos\|/,0) to project out 

components of J± . Then using Equation (3.41) 

(dE ldE\ 4neK^ - - ...... 

2 — + —^ h-V,cos(v|/ + Uz)5 (v-r), (3.45) 


(dE IdE^ 4neK^ - 
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f9(|) l9d)34jtei^^ c-3^- 
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I dz c df 1 y 


(3.46) 


The sum over all electrons can be replaced by means of an average current density 
p times the average phase, (recall \\f + k^z = ^ +cp) and Equations (3. 45) and (3.46) can 
be rewritten as 


(dE ldE3 27tei^p / ,,, 

I dz c dH y ' 


(3.47) 


where (...) is an average over sample electrons. Equation (3.47) can be rewritten in a 
complex form as 

f ^ ^ ^ _ 27te.^P / -ir\ 


1+12 W = 

dz cdt , 




(3.48) 


If the optical and electron pulses are long, then the spatial dependence in the z 
direction can be ignored and Equation (3.48) can be simplified to 




(3.49) 


39 



where a is the derivative of the dimensionless optical field a=\a\e'^ with respect to the 
dimensionless time variable x, and j is the dimensionless current density 

j - 8N{eKKLy p . Equation (3.50) is called the FEL wave equation and shows 

again that with zero average phase change and no bunching, there is no change in optical 
field. In order to achieve gain, electron bunching around ^ =7t is required. 

7. Optical Gain 

The FEL gain is defined as the fractional change of the optical power over one 
pass through the undulator. To achieve gain, energy must be removed from the electron 
beam and added to the optical beam. The change in energy, , of a single electron is 
given by 

AE^-Aymc^, (3.50) 

and is related to the change in phase velocity according to Equation (3.39) 
Av = 47t N Ay /y. Assuming a monoenergetic electron beam enters the undulator with 

initial phase velocity , the average energy change, AE^ , can be written as 

A^ = ^^}~^° ymc\ (3.51) 

" 4nN 

where (v) is the average electron phase velocity at the end of the undulator. 

Assuming a weak field (lfll«7t), low-gain (j <%) EEL, the average phase 
velocity can be approximated using a power series expansion in I a I, 

= +...). (3.52) 

Carrying out the calculations demonstrates that the average of the second order term 
produces the net change in the average energy [Ref. 19] 

= cos(vj)-l-F^vjsin(v„x) . (3.53) 

L 2 J 

Using energy conservation, the gain G can be computed as the change in energy 
of the electron beam in a differential volume dV divided by the initial optical beam 

energy in the same volume. The radiation energy in a volume element is 2E^dV /87t . 
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To obtain the total energy change of the electron beam, Equation (3.51) must be 
multiplied by the number of electrons p FdV , in that volume, where p is the electron 
density and F is the filling factor (the ratio of the electron beam area to the optical beam 
area). The gain in the optical beam is therefore 


(p FdE) (y me V 471A/^) (^ -V ^) 
lE^dVISa 


(3.54) 


and using Equation (3.53) becomes 

G=^[2-2cos(vj)-vj sin(v„x)], 


(3.55) 


where j — 8N{cKKL)^ pF ly^mc^ is the dimensionless current density (now including 
the filling factor). Equation (3.55) shows again that the optical gain in an EEE is a 
function of the initial phase velocity of the electrons. As illustrated in Eigure (3.4), the 
gain is zero at resonance (v^ = 0)- Gain reaches its maximum value when the phase 


velocity is = 2.6 (Eigure 3.5), in full agreement with our results in the phase space 
analysis. 

8. Efficiency 

The efficiency p of a Eree Electron Easer is defined as the fractional loss in 


energy of the electron beam during one pass through the undulator [Ref. 19] 

^ Ay _ i^y) 

^ 2 ’ 
ymc y 


(3.56) 


where (Ay) is the average energy of an electron at the end of the undulator. Combining 
Equations (3.38) with (3.56) gives 


(Av) 
r\ = --— 
4'kN 


(3.57) 


This chapter describes the basic theory behind SSEs and EEEs. Based on EEE 
theory, simulations were conducted to study the performance of the proposed high power 
EEE weapon design for shipboard application. 
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IV. SIMULATIONS OF fflGH-POWER FREE ELECTRON 

LASERS 


A. INTRODUCTION 

A high power free electron laser for ship defense, is being designed at the Naval 
Postgraduate School in collaboration with Jefferson Laboratory and Advanced Energy 
Systems. A short Rayleigh length resonator has been proposed to increase the spot size at 
the mirrors, and hence reduce mirror damage. This implies a nearly concentric cavity 
with a very small optical mode waist. A slight misalignment of the electron beam, either 
an offset or a tilt, could conceivably reduce the overlap between the electrons and the 
optical mode in the center of the undulator, resulting in less gain and efficiency. Three- 
dimensional simulations are used to study the effects of off-axis shifting and tilting the 
electron of the electron beam. A theoretical overview on characteristic cavity modes is 
given before describing the design parameters and the simulations methods. 

I. Optical Cavities 

In order to obtain a complete description of the beams produced by a laser, the 
wave equation must be solved. An analytic solution to the wave equation yields 
electromagnetic fields that are completely specified at all points in space. A cavity mode 
is a field distribution that reproduces itself in relative shape and in relative phase after a 
round trip through the system. A complete set of modes for a laser cavity can be 
generated by the requirement that these modes must match the mirror surfaces [Ref. 17]. 

In contrast to planar mirrors, spherical mirrors provide a more stable resonator 
configuration for the confinement of light, less sensitive to misalignment. A spherical 
mirror resonator is illustrated in Figure 4.1. 



Figure 4.1. Geometry of Spherical Mirror Resonator. Ri, R 2 are the Radiuses of 

Curvature of the Two Mirrors. 
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2. Gaussian Modes 

Gaussian beams are the modes of a spherical-mirror resonator. They provide 
solutions to the wave equation under the boundary conditions imposed by the spherical- 
mirror resonator. In a vacuum, the diffraction of the optical wavefront is described by the 
parabolic wave equation, which, in dimensionless form, is given by [Ref. 11] 

= (4.1) 

dx 4 


The Gaussian wave solution to (4.1) is given in non-dimensional form by [Ref. 9] 




a(r,x)=——e ' ' " '■'-‘’e" 
w(x) 

i.t:-xJ,, r^Cx-xJ 


w(x) 

(|)(r,x) = -tan “‘(^^—'—) + 


(4.2) 


4 + 


where r^=x^ + y^, w^(x ) = l + (x/ zl is the normalized beam width, 
[w = W/W^) and (l + z' /Z^^), is the beam waist radius, is the dimen¬ 

sionless field strength at the beam waist, x = z/L is the dimensionless time along the 
undulator, x^ is the position of the beam waist, z„=Z^/ L is the normali z ed Rayleigh 

length, and Z^-tz Wf / is defined as the distance from the beam waist over which the 
optical mode doubles in area. 



A 

V 


IT 


Figure 4.2. Expansion of a Gaussian Beam. 


As the beam moves away from the center (x^^,), it spreads out by diffraction. The 

Rayleigh length describes how fast the beam is diffracted. 

3. Resonator Optics with a Short Rayleigh Length 

At the MW power levels required for a weapon design, the resonator mirrors of 
the oscillator are vulnerable to damage because of the small beam size. Normally the 
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optical mode in a FEL has a diameter of a few millimeters (~ 4 mm). Assuming 25% 
output eoupling there is a stored power of 4 MW, so that the power density on the mirrors 
would be about 30 MW/em . No mirrors eould survive this intensity. It is estimated that 
eooled sapphire mirrors with a transmitting eoating ean handle intensity up to 200 
kW/em [Ref. 20]. A short Rayleigh length resonator has been proposed in order to 
increase the spot size on the mirrors and reduee the intensity. However, in that ease, the 
optical mode waist radius becomes very small (in the nearly concentric cavity) as seen in 
Figure 4.3. A slight misalignment of the eleetron beam, either an offset or a tilt, eould 
significantly reduce the overlap between the electrons and the optical mode in the center 
of the undulator, resulting in less gain and efficiency. 



The goal is to design the MW FEL to fit inside the 12mx4mx2m box, whieh 
means the resonator can only be 12 m long. A Rayleigh length of =1.8 cm increases 
the spot size in the resonator mirrors separated by 5 = 12 m to over 5 em diameter so that 
the intensity on the mirrors is reduced to only about 200 kW/cm . 

B. HIGH-POWER FEL DESIGN PARAMETERS 

The FEL design calls for a 185 Mev electron beam with a peak eurrent of 3.2 kA, 
electron bunches of 0.1 mm length and a repetition rate of 750 MHz. This corresponds to 
a dimensionless eurrent density of y = 210. The average beam eurrent is 0.8 A, so the 
electron beam carries 148 MW of average power. To obtain 1 MW of output power, the 
FEL extraetion efficiency must be at least 0.7 %. 

The linearly polarized undulator consists of N = 20 periods, eaeh of = 3 cm 
length, for a total length of L = 60 cm. The rms undulator parameter is = 2.8, yielding 
an optical wavelength of A- = 1 pm. The optical cavity is S = 12 m long, with a quality 

factor = 3, (2„ = l/a„ where is the % energy lost per pass) eorresponding to 28% 
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power transmission per pass. The Rayleigh length is =1.8 cm and the optical mode 
waist radius is only =0.1 mm. The PEL design parameters listed in Table 4.1 are 
illustrated in Figure 4.3. 


Electron beam 


Energy: 

i^E= 185 MeV 

Peak current: 

I peak = 3.2 kA 

Average current: 

lavg = 0.8 A 

Bunch charge: 

^ = 1.1 nC 

Pulse length: 

4 = 0.1 mm 

Repetition rate: 

Q. =750 MHz 

Radius: 

Ve = 0.14 mm 

Emittance: 

= 24mm-mrad 

Magnetic undulator 


Period: 

= 3 cm 

# of periods: 

N =20 

Length: 

L = 60 cm 

RMS parameter: 

K =2.8 

Optical Resonator 


Length: 

5= 12 

Quality factor: 

Qn = 2 

Wavelength: 

A, = 1 pm 

Rayleigh length: 

Zo = 1.8 cm 

Mode Waist: 

Wo = 0.1 mm 



Dimensionless 




parameters 




Current density: 

7 = 210 



Electron beam radius: 

Q 

II 

p 

Pulse length: 

II 

o 

Rayleigh length: 

Zo = 0.03 

Mode waist: 

Wo = 0.2 


Table 4.1. High-Power FEL Parameters. 


4=1 pm 

H H 

- 5'= 12 m -> 

Figure 4.4. MW Level FEL Design. 


E-beam : 185 MeV 
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C. BETATRON MOTION (LINEAR UNDULATOR) 

In Chapter III, the electron motion inside the undulator was studied and the 
electron trajectories for a helical undulator derived assuming perfect on-axis injection of 
the electrons. The linearly-polarized undulator, in our design, has a field of the form [Ref. 
11 ] 

= B {O,sm(k^z)cosh(k^y),cos(k^z)smh(k^y)), (4.3) 


where kg = ln / \ is the undulator wavenumber, and is the undulator period. The 

perfect electron trajectories are sinusoidal in the vz-plane, similar to the helical undulator, 
with no motion in the y-direction. Howew, when electrons drift away from the 
undulator axis, this undulator provides an additional long-wavelength oscillation and 
focusing in the y-direction called “betatron motion”. 

Assuming perfect injection in the v-direction, the relativistic force in Equation 
(3.19) can be solved giving the electron trajectories in a linear undulator [Ref. 11], 


Px =-V2—cos(^„z)cosh(fc^y), 

Y 

(4.4) 

cK^k 

P, = 2 " cos"(/r„z)sinh(2/c„y), 

Y 

(4.5) 

cK^k 

P, = " sm(2/c„z)cosh"(/r„y). 

Y 

(4.6) 


To remove the fast wiggling motion in the xz-plane, we average over an integral number 
of undulator periods . Averaging in the z direction yields |3^ = 0 ^ ~ 1. 


The averaged equation of motion for y then becomes 




(4.7) 


which describes a simple harmonic oscillator with the dimensionless betatron frequency 

KkL 2nNK 

=-. (4.i 

Y Y 
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For a typical FEL (Op « 27t so there is approximately one oscillation along the undulator. 

The general solution to the harmonic oscillator Equation (4.7) can be written as 

E9v 

y(x) = y„cos(a)pX)+—-sm(a)pX), (4.9) 

(Op 

where and 0^, are the initial position and angle of an electron at the undulator entrance 
(x =0), with respect to the undulator axis [Ref. 19]. The initial angle is related to the 

initial velocity by 0^ = y(0)/Lwhere y =dyIdx and X = z/L is the dimensionless time 

along the undulator. In Equation (4.9), Xo and ?x can also be included, but since the 
magnetic field is independent of x in Equation (4.3), any trajectory can be displaced in x 
without change. Equation (4.9) describes the electron betatron motion in a linearly 
polarized undulator with no light present. The light does not affect the betatron, but the 
betatron does affect the amplification of light. 

The effects of both types of errors, and 0^ can be compared using their 
resultant change in phase velocity. When an electron is injected off-axis by y^ or at an 
angle 0^, its average velocity is decreased relative to an electron injected perfectly. The 
corresponding change in the electron’s phase velocity is given by 
Av =-27tN -I-Y^^)/(1-I-K^) . A beam with a random spread will begin to 

decrease the EFL’s ability to bunch electrons, and thus decrease gain, when the average 
change in phase velocity becomes of order 7i [Ref 11]. 

D. SIMULATION METHODS 

To study this EEE, a three-dimensional simulation in x, y, and t is used. The 
undulator is oriented along the z-axis with the magnets normal to the y-axis. The 
simulation uses dimensionless coordinates. Eongitudinal lengths in z are normalized to 

the undulator length L, transverse lengths in x and y are divided by , and 

angles are divided by (X /. 
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The electrons are given an initial spread in positions fy:,y) and angles (9^,0^) 

determined by the beam emittance and focusing. They can also be given an offset in 
position or angle to study the effects of beam misalignment. As the electrons pass 
through the undulator, in addition to their fast wiggling motion in the vz-plane, they also 
undergo slower betatron oscillations in the yz-plane. Equation (4.9) gives the betatron 
motion for a single electron. For an electron beam, Equation (4.9) becomes [Ref. 4]. 

0 + A0 

y = iy„ + Ay)cos (tOp (x -x p)) + ^ sin (tOp (x -Xp )). (4.10) 

(Op 

In this equation, the beam misalignment is described by an offset y^ , a tilt is 
described by an angle 0^^, while the beam emittance is described by a random position 
shift of Ay and a random angular shift of The betatron frequency is given by 

Equation (4.8) and Xp corresponds to the position about which the beam is tilted. 

In the longitudinal direction, the electrons evolve in phase space according to the 
FEE pendulum Equation (3.37). An electron’s phase velocity is given by 

V - ^ - L[(k + )P^ -k]. If the electron is injected off-axis by a distance y and at angle 
0y, its phase velocity is reduced by Av =-(cOpy^-i-0^^) [Ref. 4]. 

The optical wavefront is described by the parabolic wave equation. 

The optical wavefront starts with an initial Gaussian profile, and evolves over many 
passes until a steady-state mode is obtained. The extraction efficiency is calculated using 
T] = Av /47tA^ from Equation (3.57), where Av is the average shift in electron phase 
velocity due to the FEE interaction. 

Figure 4.5 presents a 3-D simulation of the proposed laser in x, y and X . In this 
simulation, the electron beam enters the undulator normal to the xy-plane and no offset in 
position or angle is applied (y^^O. 0>o-O)- The upper right table presents the 
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dimensionless parameters describing the MW level design, along with the color scale for 
the intensity plots of the optical amplitude \a\. The top-left plot, |a(y,n)|, presents the 
evolution of a slice through the optical mode over n = 140 passes showing how the 
steady-state develops. The top-center plot, |a [x, y)|, presents the final optical wavefront 
at the undulator exit X = 1 showing the electron beam centered in the mode. The center 
plot, |a(y,x )|, depicts a section through the optical wavefront during the final pass. The 

bottom left plot,/(v,n), shows the development of the electron phase velocity distribution, 
and next to it is the final electron phase space plot at x =1. The development of optical 
power P{n) and gain G{n) is in the bottom-right. The extraction efficiency was found to 
be r| = 2.5% which is far beyond the needed value of 0.7%. 


FEL WAVEFRONTS 
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Figure 4.5. Three-Dimensional Simulation Results for the Proposed MW FEL with a 
Short Rayleigh Length and a Narrow Electron Beam Injected without an Offset in 
Position or Angle (y^ = 0, = O) ■ At Saturation, the Efficiency Reaches r| = 2.5%. 
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E. SIMULATION RESULTS 

Simulations were conducted to study the effects of electron beam misalignment. 
Specifically, the effects of tilting the electron beam are studied for two cases; tilting 


about the center of the undulator (xp =0.5) and tilting at the beginning of the undulator 



I. Tilting about the Center of the Undulator (x^, =0.5) 

The first set of simulations was run to investigate how sensitive the PEL operation 
is to electron beam tilt around the undulator center as illustrated in Figure 4.6. 



Figure 4.6. Electron Beam Tilt (Xp = 0.5). 


Eigure 4.7 shows the simulation results for steady-state efficiency r] versus initial 
phase velocity for four values of normalized electron beam tilt, 0^^ = 0, 2, 4, 6, about 
the center of undulator. The efficiency increases up to a peak, and then drops off sharply 
at the value of where the FEE gain falls below threshold. As the tilt angle is 

increased, the optimum value ofv^ increases. 
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Initial Phase Velocity, Vo 

Figure 4.7. Single-Pass Extraetion Effieieney, r| vs Initial Phase Veloeity v„ for 
Different Values of through the Undulator Center (Xp = 0.5). 

Eigure 4.8 summarizes the results of many simulations, showing the peak 
effieieney r| versus tilt angle 0^,^. As the tilt angle inereases, the overlap between the 

electrons and the optical mode waist in the center of the undulator decreases, resulting in 
a reduction in the EEE efficiency. The normalized tUt angle that corresponds to 1 mrad is 
indicated by a tickmark on the horizontal axis at 0^^ = 1.4. The efficiency begins to drop 

off for 0y^ > 3, but remains above the 0.7% goal for 0^^ < 8, well beyond the design 

tolerance of 20 p,rad , which corresponds to normalized angle 0^^ = 0.03. This relative 

insensitivity to beam tilt is due to the fortuitous choice of tilt about the center of the 
undulator, where the optical fields are most intense. 
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Electron Beam Tilt, 0^ 

Figure 4.8. Peak Single-Pass Extraction Efficiency, r| vs. Normalized Electron Beam 
Tilt 0^,^ through the Center of Undulator (Xp = 0.5). The Required Efficiency for 1 MW 
Goal (0.7%) is Achieved for Tilt Angles Eess than 6mrad. 

A single run output corresponding to electron beam tilt 0„o = 4 and initial phase 
velocity =10 is depicted in Eigure 4.9. The electron beam tilt causes the optical mode 
|a (y,x )| to tilt about the center of undulator (center plot). The optical wave amplitude is 

plotted in yellow at the mirrors on each side, indicating how toe power is distributed 
across the mirrors, and the final mode shape. The power P(n) and gain G(n) reach 
saturation after about fifteen passes (bottom-right plot). The steady-state efficiency has 
been decreased from 2.5% to 2.1%. The phase velocity is reduced by Av = -0^^^ = -16. 

This implies a change in optical wavelength of about AA-/A-= Av/27r A = 13% . The 
electron beam energy would be adjusted to compensate if this actually occurred. 
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Figure 4.9. Three-Dimensional Simulation Results for an Injeeted Eleetron Beam with 
an Initial Tilt 0^^ -4(xp =0.5). Efficiency Decreases to r] = 2.1%. 

2. Tilting at the Beginning of Undulator = O) 

The second set of simulations was run to investigate how sensitive the EEE 
operation is to beam tilt at the beginning of an undulator as illustrated in Eigure 4.10. 
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Figure 4.11 shows the simulation results for steady-state effieieney t] versus 
initial phase veloeity for three values of normalized eleetron beam tilt, 0^,^ =0, 1, 1.5, 
oeeurring at the beginning of the undulator. The effieieney again inereases ip to a peak, 
and then drops off sharply at the value of where the FEL gain falls below threshold. 

As the tilt angle is inereased, the optimum value of decreases. For (0^^ = 1, v„ = 9) 

and (0^„ = 1.5, =7), the optical mode and efficiency oscillate. A bar indicates the 

amplitude of oscillation. 



Figure 4.11. Single-Pass Extraction Efficiency, r| vs Initial Phase Velocity v„ for 
Three Different Values of 0^^ Occurred at the Beginning of the Undulator (xp = O) • 

Eigure 4.12 shows the results for peak efficiency r| versus tilt angle 0^„ when the 
electron beam tilts at the beginning of the undulator. As would be expected, the EEE 
operation is now more sensitive to beam tilt, dropping below the MW output goal of only 
0y^ > 1.5, but still well beyond the experimental design tolerance. 
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Figure 4.12. Peak Single-Pass Extraction Efficiency, r) vs. Normalized Electron Beam 
Tilt at the Undulator Entrance (Xp ~ O) • The Efficiency for 1 MW Goal (0.7%) is 
Achieved for Tilt Angles Eess than Imrad. 

A single run output corresponding to electron beam tilt 83 , 0 = 1 and initial phase 

velocity = 9 is shown in Eigure 4.13. The optical mode |a(y,n)|, and the power P(n) 

now oscillate over many passes between two different modes. The efficiency also 
oscillates between 1 % and 2 %. 
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FEL WAVEFRONTS 
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Figure 4.13. Three-Dimensional Simulation Results for an Injected Electron Beam with 
an Initial Tilt 0^,^ =1 (Xp = O). Efficiency Oscillates between 1% and 2%. 

F. CONCLUSIONS 

The required efficiency for the 1 MW goal, which is 0.7% for our design, is 
achieved in both cases for tilt angles much larger than the experimental design tolerance. 
The EEE operation is much more sensitive to beam tilt at the beginning of the undulator 
but still safely beyond the design tolerance of 20 prads. In both cases, the tilt angle 

increases the optimum value of phase velocity . This change increases for x p = 0.5 and 
decreases for (Xp = O) implying a change in optical wavelength. 

Simulations conducted to study the off-axis shift of the electron beam showed 
similar results. The MW goal is achieved for values of < 0.75 (= 0.3 mm), which is 
well beyond the design tolerance of 0.01 mm [Ref. 4]. 
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The proposed short Rayleigh length design, with a large spot size at the mirrors 
and a small optical mode waist, was found to be tolerant of electron beam misalignments 
and providing enough efficiency for the MW level weapon system requirements. 


Additional simulations, as described in the next chapter were conducted to 
explore the performance of step-tapered undulator for the upcoming 100 kW upgrade in 
Jefferson Laboratory as a transitional step from the kW level to high power (MW level). 
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V. SIMULATIONS OF 100 KW FEL USING A STEP-TAPERED 

UNDULATOR 

A. INTRODUCTION 

Jefferson Laboratory has constructed and operated the world’s highest average 
power FEL at 2.1 kW. A low risk 10 kW IR upgrade, as a next step to high power, is 
underway. All the hardware that is being installed for that system is believed to be 
capable of 100 kW operation with only slight modifications. Due to the energy recovery 
configuration with a recirculating electron beam, it is important to extract the maximum 
energy from the electron beam in a pass through the undulator while inducing the 
minimum amount of exhaust energy spread. A larger energy extraction reduces the 
requirement for a large recirculating current, while a smaller exhaust energy spread 
allows the intense electron beam to be recirculated without damaging components. 
Simulations were conducted to explore the use of a step-tapered undulator, which alters 
the resonance condition halfway through the undulator, in order to improve the FEL 
performance. 

I. Short Pulses 

As discussed in Chapter II, RE linear accelerators produce short electron pulses 
on the order of picoseconds. In general, the space-time structure of the laser beam 
approximately reflects the electron beam structure [Ref. 11]. Thus, an electron pulse with 
length 4, will produce an optical pulse with approximately the same length. The time for 
an optical pulse to complete one round-trip is determined by the length of the optical 
cavity and is equal to 25/c . To keep the system running, the reflected optical pulse must 
be synchronized with the sequence of incoming optical pulses. 

However, inside the undulator the relativistic electrons travel slightly slower than 
the photons. The resonance condition implies that during one pass through the undulator, 
the electrons will lag behind the light by the slippage distance NX . If the electron pulse 
length is comparable to the slippage distance or less, then it will not overlap the optical 
pulse over the entire pass. As the electrons slip back relative to the optical pulse, they 
amplify only the trailing edge of the optical pulse. Thus, the optical pulse centroid 
effectively travels slower than the speed of light. Eventually, the slower moving optical 
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pulse centroid will “walk” entirely off the electron pulse so that no coupling can occur 
and the optical pulse will decay due to the cavity losses. 


To compensate this effect, the optical pulse must arrive slightly ahead of the 
electron pulse at the beginning the undulator so it is not in exact synchronism. We define 
“desynchronism”, d =-2AS I NX , the displacement between electron and optical pulses 
at the entrance of undulator. In practice d is adjusted by actually moving one of the 
resonator mirrors a small distance in order to shorten the resonator and reduce the optical 
pulse round-trip bounce time. 

2. Tapered Undulator 

In the tapered undulator, the magnetic period or magnetic field strength changes 
from one end to the other. It is designed to improve the extraction efficiency, at the 
expense of reducing the small-signal gain [Ref. 9]. In contrast to the linear tapered 
undulators where the magnetic field changes gradually, a step-tapered undulator abruptly 
changes the value of the field halfway through the undulator (Figure 5.1). The step- 
tapered undulator is described by the modified pendulum equation [Ref. 5] 

n 


V = - 0 


dx‘ 


where 


X — 

V 2, 


0(z) = - 


A-l-1 a lcos(^ -l-(|)) 


(5.1) 


[0 z <0 

1 z>0 


is the step function, A = -{AtiNK^ l{\-^ K^y\{NK!K) is the amount of step-taper at 
X = 1 / 2 along the undulator, and AK /A' is the step change in the undulator parameter. 


A>0 


A<0 




nun 


Figure 5.1. Schematic for a Step-Taper Undulator 
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Previously published results have shown that as much as a 75% efficiency 
enhancement was observed with a negative, step-tapered undulator based on simulations 
and experimental results from CLIO and FELIX [Ref. 5]. 

B. 100 KW FEL DESIGN 

The 100 kW Jefferson Lab FEL design calls for a 210 MeV electron beam with a 
peak current of 270 A, in bunches of length 0.1 mm at a repetition rate of 750 MHz. The 
average electron beam power will then be 14 MW. For a 100 kW output power, an 
extraction efficiency of at least 0.7% is required. The 100 kW FEE parameters are listed 
in Table 5.1. 


Electron beam 


Energy: 

XE = 210MeV 

Peak current: 


Pulse length: 

4 = 0.1 mm 

Repetition rate: 

O = 750 MHz 

Magnetic undulator 


Period: 

\ = 8 cm 

# of periods: 

N =36 

Eength: 

L = 288 cm 

RMS parameter: 

K = 1.68 

Optical Resonator 


Eength: 

S = 36 

Quality factor: 

Qn = 4-2 

Wavelength: 

II 

B 



Dimensionless 




parameters 




Current density: 

J = 5 

Pulse length: 

II 

o 

Step-Taper rate 

A = 0,±7t ,±27U 

Desynchronism: 

d = 0 ^ 0.3 


Table 5.1. 100 kW FEE Parameters 

C. WEAK FIELD GAIN 

To study FEE gain, we use weak optical fields, where the dimensionless optical 
field is I a k 7t . The desynchronism between optical and electron pulses was varied from 
d =0 toO.3 . The gain results from many simulations with a step-tapered undulator and 

different values of desynchronism are plotted in Figure 5.2. The conventional non- 
tapered case (A = 0) produces the highest gain of 60%. At small and large values of 
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desynchronism, the gain decreases for all the undulators away from peak values around d 
= 0.1. Gain for small step-taper of A = ±3l is only slightly reduced from the conventional 
case. Larger step-taper of A = ±271 causes a significant reduction in gain to around 40% 
peaking at lower values of d ~ 0.05. The PEL still works beyond the losses for larger 
values of d< 0.2. 



desynchronism d 


Figure 5.2. Weak Field Gain G Versus Desynchronism for Step Taper with = 4.2. 
Gain Above Threshold for All Cases, Optimum Gain with No Taper (A =0). 

The results of this research also showed only a slight improvement in efficiency 
with a small negative taper (A = -7l). In all other cases explored, the efficiency is 
reduced compared to the conventional undulator (A = 0) [Ref. 5]. 

D. ENERGY SPREAD 

Figure 5.3 shows the energy spread Ay/y plotted versus desynchronism d for 
different values of A . We found that energy spread is slightly reduced for almost all 
values of positive and negative taper and stays well below the limit of 15% required for 
recirculation in all the undulator designs. 
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Figure 5.3. Energy Spread Ay/y Versus Desynchronism For Step Taper with Q^ = 
4.2. Induced Energy Spread does not Change Much with any Undulator Design. 

E. CONCLUSIONS 

The results of this research showed that step-taper is not always effective in 
increasing the FEE performance. An FEE must be far into the strong-field saturation 
before tapering can significantly extend the saturation limit. Only slight increase in 
efficiency over the conventional undulator was observed for small negative taper. The 
induced energy spread did not change significantly for any undulator design and was 
sufficiently small for recirculation. However, in all cases the gain was found to be above 
the threshold of 21% required for the 100 kW design. 


Another important issue for operating a High Power Easer is the atmospheric 
propagation of the laser beam. The atmospheric effects can put more strict limits to the 
power requirements of the DEW systems, as described in the next chapter. 
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VI. POWER SPECTRUM OE THE EEL 


Earlier studies have shown that in order to destroy an oncoming missile, the laser 
power required on the missile surface is about 12 kW/cm^ over a radius spot of Wo = 5cm 
for a two-second duration [Ref. 18]. These numbers determine the minimum output 
power (~ 1 MW) required in order to achieve an effective lethal EEL weapon system 
(ignoring the propagation losses). However, another issue of interest is how this power is 
to be distributed over the frequency spectrum. The extremely short (picosecond) pulse 
length of the EEL optical beam, a result of the short electron pulses, produces frequency 
components around the fundamental (X=lpm) with significant amounts of energy. 
This may affect the performance of the EEL as a weapon in many ways. The 
concentration of energy at certain frequencies can dramatically increase the atmospheric 
absorption and scattering, and remove a significant amount of power from the beam. On 
the other hand, microwave power may result in additional damage to the target, 
destroying the electronics and neutralizing the guidance circuits of the missile. Einally, 
the effects on the visibility of the laser beam and eye-safety issues may arise. 

A. FOURIER TRANSFORM OF A LASER PULSE TRAIN 

Consider the output of an EEL weapon consisting of a number of optical pulses 1 
psec long separated by T = I/O where O = 750 MHz is the pulse repetition frequency 
as illustrated in Eigure 6.1. The optical wavelength of the pulses is 1 pm and the total 

average power 1 MW. Each pulse contains several hundred wavelengths of light. 



T = 1.3 ns 


Eigure 6.1. EEL Output Pulse Train. 
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The mathematic representation of the pulse train is 


- nT)cos{V)J +(|) (n)), 

n 


( 6 . 1 ) 


where a{t) is the amplitude of the electric field of the pulse with peak amplitude a, fl (t) 
is the rectangular amplitude window 


n(o = 


1 


[0 


T 1 
Ul< — 

2 k 

otherwise 


( 6 . 2 ) 


n is the number of pulses, n) is the random phase of the optical wave in each pulse, 
0)^ = 27 c/„ and f^=c/X is the frequency of light. Since the phase ^(n) does not 
contribute to the energy, and hence, to the time-average power of a(t ), it will be ignored 
for convenience sake in the calculations and a pulse train of n identical pulses considered. 
Thus 


a(t) = aY,U(t -nT)cos( 0 )J ), 

n 


(6.3) 


a periodic signal with period T. 

Since a 2 second burst contains 1.5 billion pulses for our problem, an infinite 
number of identical pulses will be considered. 

In order to discover how the power is distributed in the frequency domain, the 
Fourier series coefficients of our periodic signal a(t) needs to be computed. For 

convenience, only one pulse a(t) is considered and the Fourier transform A(/) of a 
single pulse computed. Then 

A(f) = F[fl( f)] = F[fl n()cos((o„0] • 

Taking the Fourier transform obtains 

Mf) = ^aT [sinc(/ - fo)T^ + sinc(/ -H fo)T ]. 

Based on the fundamental Fourier theory, it is known that any aperiodic signal 
can be viewed as a periodic signal with an infinite period. More precisely, as the period 
becomes infinity, the frequency components form a continuum and the Fourier series sum 
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becomes an integral. The relationship between the Fourier Transform of a(t) and the 


Fourier coefficients of a{t) is given by [Ref. 21] 

(6.4) 

or 

a, = ^[sinc(/ - /o)r„+sinc(/ + . (6.5) 

From Parseval’s relation for continuous-time periodic signals, the average power 
of the periodic signal is [Ref. 21] 

= f, 

k=-oo 

SO that lOC^ f is the average power in the km harmonic component of a(t). 

The above equation can also be written as 

+“ fjT 

K. = E ' - /•)r +sinc(/ + m F|^_^ , (6.6) 

or 

+“ a T 

P.. = 

k=0 

since the power exists only in positive frequencies. 

Moreover, the time-average power Pqv of a{t) during the time interval T is given 
by 

P,.=Y' 

where now 

E. = K'; 'a cos* (271 fj)dl , (6.9) 

and the sum in Equation (6.7) evaluates to 
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(6.10) 


P 

Substituting the numbers in the above equation, the amplitude a of each pulse can 
be calculated. Equation (6.7) with a simple code can now be implemented and the power 
spectrum of the laser beam studied. 

B. POWER DISTRIBUTION 

The power distribution over the frequency spectrum is depicted in Figure 6.2 (a). 
Most of the power is concentrated around fo = 3x10*"^ MHz which corresponds to 
A, = 1 pm, while the average power substantially cfecreases when moving to higher or 
lower harmonic components. The main lobe (Figure 6.2 (b)) contains 90% of the total 
power of the laser beam, P^i ~ 0.9 MW. 




Frequency x (Hz) 

(b) 

Figure 6.2. Power Spectrum around 1 pm (a). Main Fobe (b). 
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Thus, a huge amount of power (= 10^ W) is contained in other frequency 
components as illustrated in Figure 6.3. 

19 17 

The frequency region between/= 10 Hz and/= 10 Hz contains almost 99.99% 
of the total power of the beam while the power in the harmonics outside this region is 
negligible. 



Frequency x (Hz) 

Figure 6.3. Normalized Power Spectral Density of the FEL Beam. 

Specific frequency ranges of interest can now be further studied. 

1. Microwave Regime 

The microwave region is very important as a significant amount of power in 
microwaves could cause a catastrophic failure to control and navigation electronics of the 
target missile. However, the power contained in the microwave region (A- = 1 mm to A 
= 1 m) is found to be too little to result in a kill mechanism, Pmic = 0.085 Watts (Figure 
6.4). Thus, the FEL weapon is not expected to seriously affect and disorient an incoming 
cruise missile through the microwave region of its spectrum. 
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Figure 6.4. Power Spectra in Microwave Regime. 

2. Visible Spectrum 

Contrary to the microwave regime, the visible region is located very close to the 
near infrared, where the fundamental (1 pm) is located, and thus, a significant amount of 
power is contained in this region. Figure 6.5 depicts the power spectral density in the 
visible region =400nmto^ =780nm). The total power contained in the visible 
region is found to be Pyis ~ 500 W for the square optical pulse. 


t 10 



Figure 6.5. Power Spectral Density in the Visible Region. 
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Now, it can further estimated if the PEL beam will be visible to the naked eye. In 
the configuration of Figure 6.6, a 1 MW FEL burst against an oncoming missile 5 km 
away from the beam director is considered. 



Figure 6.6. Scattered Power in the Visible Spectra. 


First, it is necessary to compute the total visible power scattered through the 
atmosphere along the beam path. In the optimum case of no clouds or fog, the scattering 
rate in the visible region is Gsc = 0.1/ km (Figure 6.7). The scattered power is given by 
[Ref. 22] 

Psc = Pvis(l-e~‘’-^), ( 6 . 11 ) 

where Psc is the scattered power in the visible region and Pvts = 500W. For z = 5 km, Psc 
= 200W is obtained. 


Assuming that the light scatters uniformly in all directions, the intensity of the 
scattered light is given by 



Sfo 2k rz 


( 6 . 12 ) 


where r is the average radius of the beam ( r = 10 cm) and St the total external surface 

9 9 

area of the beam. From Equation 6.12, Pc = 0.06 W/m or 60 nW/mm . For the human 
eye, the radius of the pupil is 2 mm and the pupil area 3 mm [Ref. 23]. Assuming that 
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the observer is stantioned right next to the beam path, the fraction of light entering the 

'y 

eye is then Peye = Pc mm =180 nW. 

If the divergence of scattered light as well as the propagation losses between the 
PEL beam and an observer are now considered, it is thus observed that the power level 
entering the eye is even lower and unlikely to be detected at any reasonable distance 
(even several meters). 



Figure 6.7. Aerosol Scattering. (From [Ref. 1]) 


C. ATMOSPHERIC EFFECTS 

Many of the challenges in developing a HEL weapon system that will operate in a 
maritime environment are related to atmospheric effects and propagation. Of primary 
concern to the Navy are the effects of high aerosol scattering on the propagation of the 
HEL beam from the source to the target. Natural aerosols including fog, precipitation, 
and clouds must be considered in any planned use of laser weapons. On the battlefield, 
airborne dust thrown up by exploding ordnance and smoke from burning equipment, may 
also occur. 


On the other hand, absorption by atmospheric gases and aerosols can generate 
atmospheric density gradients that cause the laser beam to spread out. This effect is 
known as thermal blooming. 

It is generally understood that scattering effects become more significant as the 
laser wavelength shortens (Figure 6.7). However, scattering is a linear process and 
depends on the index of refraction and the mean particle size. In contrast, absorption 
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leads to non-linear losses for high-energy lasers, where an inerease in power from the 
laser results in an inereasing pereentage of energy loss. Thus, an atmospherie window 
with low absorption is important for minimizing thermal blooming. 

At the moment. A, = 1 pm seems to be the best compromise wavelength between 
scattering and blooming and has been chosen as the operating PEL wavelength. 
Absorption and scattering are the primary causes for loss of energy in the atmosphere, but 
they must be studied in detail to understand their effects on weapon system performance. 

1. Power Removal (Pulse Distortion) 

The energy losses due to aerosols are described by where Ue is called the 
“extinction coefficient” and represents the aerosol attenuation due to both absorption and 
scattering as a function of wavelength. In reality, the maritime conditions at a particular 
geographic location may vary on different days. The extinction coefficient has been 
described by computer modeling codes such as FASCODE, EOWTRAN, MODRAN, 
and MOEEY. 

All of these codes allow the user to insert weather conditions and select from 
different environments. It is obvious that different conditbns can lead to different 
results. For the purpose of this work, a specific maritime aerosol model for a moderately 
clear atmosphere will be considered (Figure 6.8). 



Frequency x (Hz) 

Maritime Aerosol Model (23 km Visibility). (After [Ref. 24]) 
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Figure 6.8. 





It was demonstrated previously that the short pulses of the FEL beam lead to a 
broadband spectrum with significant power in harmonics far from the fundamental (A- = 
1 pm), where each one has a different value of extinction. 

Since the power contained in each wavelength is known, the aerosol attenuation 
coefficients can be applied to the FEL beam power spectra and the total power P,nis that 
can be delivered from a laser to the target computed 

, (6.13) 

X 

where Pship (A ) is the laser beam output power corresponding to each wavelength. 

The amount of power that is delivered through the atmosphere over a distance z = 
5 km using the specific model is found to be P,„is = 0.47 MW. Thus, even though 1 MW 
is enough to destroy an incoming missile, atmospheric effects result in significant 
reduction in on-target irradiance so that lethal effects may not be obtained (Figure 6.9). 

The total power on the target could be actually increased if the FEL operated at 
shorter wavelengths. As can be seen in Figure 6.7, by shifting to bnger wavelengths 
(lower frequencies), the main lobe, which contains the 90% of the total power, can be 
moved to smaller values of extinction, since this is basically determined by the scattering 
coefficient. However, when moving to regions with high atmospheric absorption the 
problem of thermal blooming arises. 



Figure 6.9. Power Removal. 
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2. Thermal Blooming 

During an engagement, the path of the laser beam column may include a region of 
stagnate air. The existence and the length of a stagnation z)ne cannot be predicted 
precisely, but statistically it is expected to be between 1 m - 100 m long [Ref. 24]. The 
thermal blooming effect is created by molecular absorption when the laser beam 
propagates through such a stagnation zone in the atmosphere. The absorbed power in this 
zone leads to temperature rise and, thus, to a change in the index of refraction. 
Accordingly, a difference in the index of refraction changes the speed of light and, thus, 
the optical path length across the beam wavefront. In other words, light travels slower 
due to low temperatures, higher air density on the outer portion of the laser beam and 
faster due to higher temperatures, lower air density, inside the laser beam. This effect, 
which is analogous to a negative lens, spreads the laser beam and can cause catastrophic 
phenomena during high-energy laser beam propagation. 

Thermal blooming is a non-linear effect. Since the temperature rise is 
proportional to the absorbed power, it can be seen that more power leads to more 
blooming. Minimizing the absorption is important to the performance of the laser as a 
weapon. 

The amount of power absorbed from the atmosphere can be computed by 
knowing the power spectrum of the laser beam. The limiting consequences of thermal 
blooming may then be evaluated for an PEL weapon system [Ref. 25]. 

In closing this section, it must be mentioned that absorption and scattering are 
generally incorporated in analysis using statistical averages. In reality, maritime 
conditions, at a particular geographic location and at a particular time, may favor 
different wavelengths on different days. The solution to this problem can be the use of an 
PEL, which even at high power levels, could be pre-designed to operate at different 
wavelengths. This option which is not possible with solid-state or chemical lasers. 

D. GAUSSIAN SHAPE PULSES 

It has been demonstrated how the power of an PEL beam is distributed throughout 
the frequency spectrum as well as how the power spectrum can provide important 
information about the atmospheric propagation effects. However, in our calculations, a 
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laser beam consisting of rectangular micropulses has been considered, since the 

amplitude a is a positive constant. In reality, Gaussian-shaped, or smooth electron 
pulses from the accelerator will produce Gaussian-shaped, or smooth, or similarly optical 
pulses. 

At this point, it is instructive to consider the effect of the shape of the micropulse 
on the power spectrum of the laser beam. The discontinuities that appear in the 
rectangular pulse shape result in high sidelobes in the frequency domain. In contrast, 
pulses that do not contain abrupt discontinuities have correspondingly low sidelobes in 
their frequency spectra. 

In order to observe the differences, the power spectral density of two micropulses 
with different shapes will be computed. To represent a Gaussian shape pulse, which 
actually extends from - 0 ° to - 1 -°°, a very common pulse shape in digital filtering design, 
the Hanning pulse, can be used instead as can be seen in Figure 6.10 [Ref. 26]. 


xio"^ Single laser pulse 



Time X (sec) ^-12 

' ' X10 


Figure 6.10. A CW Pulse with (a) Rectangular Envelope (in Blue), (b) Using a 

Hanning Pulse (in Red). 
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The “bell-shape” pulse is twice as long as the rectangular one and contains no 
discontinuities while falling off smoothly. The power spectral densities of both 
micropulses are illustrated in Figure 6.11. 


0 



Frequency x (Hz) 

Figure 6.11. Normalized Power Spectra Density of a micropulse (a) Rectangular 
Envelope (in Blue), (b) Hanning Pulse (in Red). 

It can be seen that using the Hanning pulse shape significantly lowers sidelobes 
compared to the rectangular window, while more power is concentrated into the main 
lobe or around the fundamental frequency. This means that, in reality, there will be less 
power in the harmonics far from the fundamental than computed earlier considering the 
rectangular pulses. The concentration of the beam energy around the selected frequency 
is obviously desirable and expected to ameliorate the negative effects of atmospheric 
absorption and aerosol scattering. Rectangular envelope pulses correspond to the worst 
case as far as the broadening of the spectrum is concerned, and can just provide 
indicating values of the power spectrum of the FEE beam. 

E. CONCLUSIONS 

This chapter shows how the power of the FEE is distributed over the frequency 
spectrum. The atmospheric absorption and scattering for the specific maritime aerosol 
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model showed a power removal of 53% over a distance of 5 km. However different laser 
pulse shapes can give different results since as the pulses gets smoother more power is 
concentrated on the main lobe around the fundamental frequency. The next Chapter 
outlines the conclusions of the whole research. 
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vn. CONCLUSIONS 


High-energy lasers provide the Navy with unique opportunities to augment or 
improve operational capabilities in a variety of mission areas. The most likely candidates 
for the development of a Navy high-energy weapon system are the high power Free- 
Electron Lasers (FELs) and Solid-State Lasers (SSLs) for three reasons: 

• These are both all-electric lasers, and there are inherent dangers associated 
with the storage of the chemicals needed for chemical lasers on board 
ship; 

• The need to develop a source at the appropriate wavelength for maritime 
propagation, and present-day chemical lasers are at the wrong wavelength; 

• The Navy’s recent decision early in 2000 to adopt an electric drive 
propulsion system for DD-21 and follow-on ships is consistent with these 
two all-electric lasers. 

In general, FELs can produce powerful, highly coherent, optical beams, but also 
place extreme demands on the resonator mirrors. SSLs, typically, produce a larger beam 
spot on the mirrors, requiring optics with a lower damage threshold, but they generally do 
not achieve the desirable beam quality. Efficient heat removal from the laser media is 
needed in high-power SSL designs, while high-power FELs require a high-brightness 
injector and novel optical cavity designs. However, with adequate research development 
efforts, it appears that both technologies have the potential to meet the needs of a 
weapon-class laser. 

Fiber optic lasers and the Heat Capacity Laser project are examples of the 
readiness of solid state laser technology for rapid maturation. On the other hand, the 
unique advantages of FELs, such as wavelength selectability and tunability, high beam 
quality, the potential to achieve high peak power without distorting the lasing media and 
long run-time, could favor FELs for weapon-class lasers in the near future. 

At the Naval Postgraduate School, a high power free electron laser for ship self- 
defense is being designed with a short Rayleigh length, which results in large spot sizes at 
the mirrors to avoid mirror damage. Three-dimensional simulations used to study the 
performance of the proposed design demonstrated that it is tolerant of electron beam 
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misalignments, and provides high enough extraction efficiency for the MW level weapon 
system requirements. 

The required 0.7% efficiency for the 1 MW goal is achieved in simulations for tilt 
angles much larger than the experimental design tolerance. It is found that PEL operation 
is much more sensitive to beam tilt at the beginning of the undulator than tilting around 
the middle, but it is still safely beyond the design tolerance. 

Simulations were also conducted for the 100 kW PEL design as a transitional step 
to high power. The results of this research showed that weak field gain exceeds the 
required threshold of 21%, while the strong-field energy spread did not change 
significantly over any undulator design. 

Pinally, the spectrum of the proposed PEL output was studied. The power is 
distributed over a wide frequency range. It was found that 90% of the total power is 
concentrated in the main lobe (around 1 pm wavelength), while power decays sharply 
for higher or lower harmonics. Atmospheric absorption and scattering were also studied 
for a specific maritime aerosol model showing power removal of 53% up to a distance of 
5 km. However, the effects of Gaussian-shaped pulses showed that there will be less 
power in the harmonics far from the fundamental than computed considering the 
rectangular pulses. The concentration of the beam energy around the selected frequency 
is obviously desirable and expected to ameliorate the negative effects of atmospheric 
absorption and aerosol scattering. 
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